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SUMMARY 

A computer program MCFLARE that uses Monte Carlo methods to simulate solar 
flare occurrences during an interplanetary space voyage is described. The total biolog- 
ical dose inside a shielded crew compartment due to the flares encountered during the 
voyage is determined. The computer program evaluates the doses obtained on a large 
number of trips having identical trajectories. From these results, a dose D having a 
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probability p of not being exceeded during the voyage can be determined as a function of 
p for any shield material configuration. 

The user of the code selects any number of solar flares considered to be represent- 
ative of the ones that will occur during future solar active periods (these flares are gen- 
erally selected from the flares that occurred during the last solar active period 1956 to 
1962). The flares are assumed to occur randomly during these periods. The dose at a 

O 

distance of 1 astronomical unit (1. 496x10 km) from the Sun from each of these flares 
behind any shield configuration investigated is input to the MCFLARE code. The code 
accounts for the dependence of the dose received from a flare on the distance from the 
Sun according to a (l/r) a variation, where r is the distance from the Sun and the ex- 
ponent a can be assigned any integral value including zero. From trajectory parame- 
ters, which are input to the computer program, the distance from the Sun as a function 
of time during the trip is calculated. 

To illustrate the use of the code, a trip to Mars and return is calculated, and esti- 
mated doses behind several thicknesses of aluminum shield and water shield are pre- 
sented. A FORTRAN IV listing, data input instructions, and sample output are given. 



INTRODUCTION 


The protons emitted by solar flares are the most hazardous source of space radia- 
tion encountered on an interplanetary space voyage which rapidly traverses the Van Allen 
belts. To define shielding requirements for biological protection from solar flare pro- 
tons encountered in manned interplanetary missions, a computer program MCFLARE 
has been prepared which uses Monte Carlo methods to simulate solar flare occurrences 
during the trip and records the total associated biological doses for the trip behind vari- 
ous shield material configurations. 

A set of flare events, considered to be representative of those that will occur during 
future solar active periods, is required. Generally, this set is selected from the flares 
that occurred during the last solar active period from 1956 to 1962. The flares are as- 
sumed to occur randomly during the trip (which is assumed to take place during the solar 
active period). The doses behind the various shield configurations investigated have to 
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be evaluated for each flare, at a distance of 1 astronomical unit (1. 496x10 km) from the 
Sun, by using another computer code such as one of those described in references 1 or 2. 
(Reference 3 presents doses from solar flare proton spectra behind two shield materials 
as calculated by using reference 1. ) These dose values are input to the computer pro- 
gram MCFLARE. This program provides for including the dependence of the magnitude 
of the dose received from the flare on the distance r from the Sun. This dependence is 
taken to be of the form (r^r)®, where r Q is the distance between the Sun and the Earth, 
and the exponent a can be any arbitrarily selected integer including zero. 

The Monte Carlo procedure (analogous to particle transport Monte Carlo) is used to 
determine when flares are encountered on the trip, which of the input flares occur, and 
the total dose obtained from all flares encountered during the trip behind each shield 
configuration. The computer program tallies the doses obtained for a large number of 
trips having identical trajectories and from this, the distribution of doses incurred on 
the mission for each shield configuration is obtained. Each dose distribution determines 
a dose D having the probability p of not being exceeded as a function of p. 


METHOD OF ANALYSIS 

Any N flare events can be selected as representative of the flare activity that will 
occur during the space voyage. The doses from each flare at a distance of 1 astro- 
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nomical unit (1. 496x10 km) from the Sun behind one or more shield material config- 
urations investigated are input to the computer. Each flare is assumed to occur ran- 
domly on the average of once during the solar active period T re j, which is of the order 
of 2000 days or 5^ years. 


2 



Monte Carlo Method to Determine Occurrences of Flares During Trip 


Inasmuch as flares are assumed to occur randomly in the solar active period, the 
probability of any flare occurring in the time interval dt is (N/T ,;)dt = ji dt. Then 

iif x GX 

the probability of no flare occurring for time t is e p and the probability of a flare 
occurring in dt at time t is e“^V dt. 

In the Monte Carlo method, this probability of flare occurrence in the time interval 
dt at time t is equated to the probability of selecting a random number in an interval 
d£ about £. If a random number set uniformly distributed in the interval (0, 1) is se- 
lected, then the probability of selecting a random number in d£ about £ is d£. Equat- 
ing these two probabilities, one gets 


e~ M V dt = d£ 

f l e" Mt V dt' = d|» 

Jo Jo 

t - - - ln(l - V) 

Inasmuch as both | and 1 - £ are uniformly distributed in the interval (0, 1) an 
equivalent expression for t which results in the same exponential distribution is 

t = - - In 4 (1) 

M 

By selecting a random number from a uniformly distributed set in the interval (0, 1) the 
time between flare events is obtained from equation (1). 

When an event has occurred, a particular flare of one of the N is also selected by 

IV 

random numbers. The i in flare occurs when a new random number occurs in the in- 
terval 


3— 1 < 4 i = 1,2,. . ., N (2) 

N N 

Dose values input to the code are calculated by using the methods of reference 1, or 
2, for each representative flare using the solar flare proton intensity and spectral shape 
as observed at the Earth. The code permits an effect of the distance from the Sun on the 
dose (or intensity) received from a flare of the form {r^r) a where r Q is the distance 
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from the Sun to the Earth, r is the distance of the spacecraft from the Sun during the 
flare, and a is an exponent that is assigned any interger value including zero. 

The code evaluates r for any time during the trip by calculating the mission trajec- 
tory from input parameters describing the two transfer ellipses and the total trip time. 
The first ellipse is the trajectory from the Earth to the planet and the second describes 
the return trip. Details of determining position as a function of time are presented in 
appendix A. 

The calculation proceeds as follows: Selecting a random number and using equa- 
tion (1) determines the time in the mission when a flare occurs. Selecting another ran- 
dom number and using equation (2) determines the particular flare which occurs. Know- 
ing the time in the mission (and, hence, distance from the Sun) and the particular flare, 
the dose received through each shield configuration is tallied. The duration of the flare 
event is then added to the time in the mission when the flare occurred to give the time 
when the event is over. From a new random number and equation (1), the time elapsed 
until the next flare occurs is calculated, and so on. This procedure is repeated until the 
total trip time has elapsed. The total dose from all flare events encountered on the trip 
is determined for each of the shields considered. 


Selection of Dp 

This procedure is repeated for a large number M of trips of identical trajectory. 
The code divides the dose range into dose intervals that are 1 dose unit wide and tabu- 
lates the number of trips A(D) during which doses that lie within each dose interval are 
encountered. The cumulative fraction of the total trips that encounter doses less than 
the upper dose bound of each interval is also tabulated (this is an estimate of the proba- 
bility p that a dose equal to the upper dose bound will not be exceeded). From this tab- 
ulation, for any specified p, a dose D p which has the probability p of not being ex- 
ceeded on the trip can be selected for any shield. The standard deviation of this value 
from the true value is shown in appendix B to be 



i fW 1 

y Mf 2 (D p ) 


(B6) 


where M is the number of trips in the group investigated, and f(D p ) is the dose proba- 
bility density at D p (which is calculated from the tabulated output). If A(D p ) is the 
number of trips contained in the interval in which D p occurs, f(D p ) is approximately 
given by 
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«y = 


Ady 

M 


If A(D) is a widely fluctuating value in the vicinity of D p , averaged values of A(D) can 
be hand calculated, and f(D p ) is obtained from the averaged value of A(D p ). Note that 
for a given p, because f(D p ) is independent of the number of trips, cr^ varies in- 
versely as >/m. ^ 

To reduce the uncertainty associated with D p for a given number of trips, the code 
is written so that the trips constitute a stratified sample for the time-to-first-flare se- 
lected in accordance with the exponential distribution exp[-(Nt/T re £)]. This technique 
is described in appendix C„ 


ILLUSTRATIVE EXAMPLE 

An interplanetary trip to Mars and return has been evaluated in order to illustrate 
details of the method and to present some representative shield requirements for such 
a mission. Water and aluminum were the shield materials considered. 


Selection of Flare Events 

The flare occurrences during the last solar active period 1956 to 1962 are assumed 
to be representative of those that will occur during a future active period. Webber in 
reference 4 has compiled a record of these flares, their time-integrated proton inten- 
sities, and spectra. From these records, the 20 largest flares, based on proton inten- 
sity, were selected (the effect of neglecting the rest of the flares is small for shields 

considered here). These flares are listed in table I which presents the integrated flux 

-12 

of protons having energies greater than 30 MeV (4. 8x10 J) and some constants A 
and P Q associated with each flare. The spectral shape of each flare is assumed to 
vary as 


N>E = Aexp^- V E 2 y 87 ig^ 

where N > E is the total number of protons per square centimeter with energies 
greater than E. The values of P Q were obtained from Webber’s compilation, and the 
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TABLE I. - FLARES (FROM SOLAR ACTIVE PERIOD 1956 TO 1962) 
SELECTED AS REPRESENTATIVE OF THOSE THAT WILL 
OCCUR DURING FUTURE ACTIVE PERIODS 


[Integral flux above 30 MeV (4. 8x10" 12 J), and spectral constants A 
and P Q for each flare are tabulated. ] 


Flare 

N > 30 MeV 

1 p_ 


A, 

date 

(4. 8X10“ 12 J), 
protons/ cm 2 

O 


protons/ cm^ 

2-23-56 

1.0X10 9 

195 

3.4X10 9 


1-20-57 

2X10 8 

61 

1. OxlO 9 


8-29-57 

1.2X10 8 

56 

8.6X10 9 


10-20-57 

5xl0 7 

127 

3. 3xl0 8 


3-23-58 

2.5X10 8 

64 

1. lxlO 9 


7-7-58 

2.5X10 8 

62 

1.2X10 9 


8-16-58 

4xl0 7 

64 

Q > 

1. 7xl0 y 


8-22-58 

7X10 7 

56 

5. OxlO 9 

f Clustered flare event 

8-26-58 

1. 3xl0 8 

51 

1.2X10 10 , 


5-10-59 

9.6X10 8 

84 

1. 7xl0 10 


6-13-59 

8. 5X10 7 

CO 

a 

1.2X10 10 


7-10-59 

1.0X10 9 

104 

1. OxlO 10 ^ 

1 

7-14-59 

1. 3xl0 9 

80 

2. 6xl0 10 

v Clustered flare event 

7-16-59 

9. 1X10 8 

105 

8. 9xl0 9 j 

1 

9-3-60 

3. 5xl0 7 

127 

2.3X10 8 


11-12-60 

1.3X10 9 

124 

9 'l 

8, 9xlO y 


11-15-60 

7.2X10 8 

114 

5.9X10 9 

^ Clustered flare event 

11-20-60 

4. 5xl0 7 

118 

3.4X10 9 > 


7-12-61 

ft 

o o 

CO -J 

56 

2.9X10 9 ] 
3. lxlO 9 J 

► Clustered flare event 

7-18-61 | 

102 



a Estimated. 


values of A were calculated to be consistent with flux values of N > 30 MeV 
(4. 8X10“ 12 J). 

It has been observed that some of the flares tend to occur in clusters. Table I indi- 
cates that such clusters occurred in August 1958, July 1959, November 1960, and July 
1961. These four clusters were selected as representative of the clustered flare events 
that may occur and were assumed to be events of 12-, 8-, 10-, and 8-day durations, 
respectively. The other nine single flares were each assumed to have a duration of 
2 days. During future solar active periods, each of the four clustered events and each 
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of the nine single flare events are assumed to occur randomly with a frequency of once 
during the active period (taken to be 2000 days). 


Doses For Selected Flare Events 


Table n presents the doses obtained behind various shield thicknesses from each of 
these 13 events. These doses were calculated by using the Lewis Proton Shielding Code, 
described in reference 1, and are representative of doses received at the center of a 
spherical crew compartment having the given shield thicknesses when the vehicle is lo- 


TABLE H. - RAD (cj/kg) AND REM DOSES FROM SELECTED FLARE EVEN IS EErtlND VARIOUS THICKNESSES OF WATER 
AND ALUMINUM SHIELD AT 1 ASTRONOMICAL UNIT (1. 496xI0 8 km) FROM SUN 
(a) Rad (cj/kg) dose behind water shield (b) Rem dose behind water shield 



Shield thickness, 


Flare date 



g/cm 2 




10 

15 

20 

30 

40 




Dose, 





rad (or cj/kg) 



38.32 

25.00 

18. 07 

10. 66 

6.95 

2-23-56 

.90 

.29 

. 12 

.03 

.01 

1-20-57 

.35 

. 10 

.04 

.01 

.00 

8-29-57 

1.08 

.58 

.35 

. 16 

.09 

10-20-57 

.94 

.31 

. 13 

.03 

.01 

3-23-58 

1.06 

.34 

. 14 

.04 

.01 

7-7-58 

.85 

.26 

. 10 

.03 

.01 

Aug. 1958 cluster 

9.93 

4. 15 

2.09 

.71 

.31 

5-10-59 

.49 

. 15 

.06 

.01 

.01 

6-13-59 

44.01 

20.39 

11. 19 

4.37 

2.07 

July 1959 cluster 

.76 

.41 

.25 

. 11 

.06 

9-3-60 

44.32 

23.31 

14. 13 

6.35 

3.34 

Nov. 1960 cluster 

4.88 

2.31 

1.28 

. 50 

.24 

July 1961 cluster 


(c) Rad (cj/kg) dose behind aluminum shield 




Shield thickness, 



Flare date 




g/cm 2 





10 

| 15 

| 20 

j 30 

| 40 

50 

| 60 





Dose, 







rad 

(or cj/kg) 




48. 11 

32.23 

23. 84 

14.55 

9.82 

7.02 

5.24 

2-23-56 

1.71 

.60 

.27 

.09 

.04 

.03 

.02 

1-20-57 

.71 

.23 

. 10 

.03 

.02 

.01 

.01 

8-29-57 

1.52 

.85 

. 54 

.26 

. 15 

.09 

.06 

10-20-57 

1.78 

.63 

.28 

.09 

.04 

.03 

.02 

3-23-58 

2.01 

.71 

.32 

. 10 

.05 

.03 

.02 

7-7-58 

1.69 

. 55 

.24 

.07 

.04 

.02 

.02 

Aug. 1958 cluster 

16. 15 

7. 19 

3.84 

1.45 

.69 

.39 

.25 

5-10-59 

.99 

.32 

. 14 

.04 

.02 

.01 

.01 

6-13-59 

67. 55 

32.90 

18.94 

7.97 

4.10 

2.41 

1.57 

July 1959 cluster 

1.07 

.59 

.38 

. 18 

. 10 

.06 

.04 

9-3-60 

62.93 

34. 57 

21.79 

10. 42 

5. 83 

3.62 

2.40 

Nov. 1960 cluster 

7.40 

3.67 

2. 14 

.92 

.48 

.28 

. 18 

July 1961 cluster) 



Shield thickness, 


Flare date 



g/cm 2 




10 

15 

20 

30 

40 




Dose, 






rem 




42. 10 

27.76 

20. 79 

13. 03 

9.03 

2-23-56 

1.07 

.38 

. 19 

.07 

.04 

1-20-57 

.43 

. 15 

.07 

.03 

.01 

8-29-57 

1.20 

.65 

.42 

.21 

. 12 

10-20-57 

1. 12 

.40 

. 19 

.07 

.04 

3-23-58 

1.26 

.45 

.22 

.08 

.04 

7-7-58 

1.03 

.35 

. 17 

.07 

.04 

Aug. 1958 cluster 

11.30 

4.93 

2.67 

1.09 

. 57 

5-10-59 

.60 

.20 

. 10 

.04 

.02 

6-13-59 

49.59 

23.68 

13.76 

6. 11 

3.36 

July 1959 cluster 

.84 

.46 

.29 

. 15 

.09 

9-3-60 

49.34 

26.49 

16.79 

8. 30 

4.85 

Nov. 1960 cluster 

5.48 

2.67 

1. 57 

.71 

.39 

July 1961 cluster 


(d) Rem dose behind aluminum shield 




Shield thickness, 



Flare date 




g/cm 2 





10 

| 15 

| 20 

30 

40 

50 

60 





Dose, 








rem 





53.93 

37.20 

28.77 

19.07 

14. 16 

11. 15 

9.26 

2,23-56 

2. 11 

.86 

.48 

.24 

. 16 

. 12 

. 10 

1-20-57 

.90 

.35 

. 19 

. 10 

.07 

.06 

.05 

8-29-57 

1.73 

1.00 

.68 

.38 

.25 

. 18 

. 14 

10-20-57 

2.20 

.89 

. 50 

.25 

. 17 

.12 

. 11 

3-23-58 

2.48 

1.01 

. 56 

.28 

. 19 

. 14 

. 12 

j 7-7-58 

2.12 

.83 

.45 

.24 

. 17 

. 13 

. 11 

Aug. 1958 cluster 

18.98 

9.06 

5.36 

2.59 

1.63 

1. 18 

.93 

5-10-59 

1.25 

.49 

.26 

. 14 

. 10 

.08 

.06 

6-13-59 

78.42 

40.39 

25. 19 

12.79 

8. 12 

5. 86 

4. 59 

July 1959 cluster 

1.22 

.71 

.48 

.26 

. 17 

.13 

! . 10 

9-3-60 

71.89 

41. 17 

1 27. 59 

15.09 

9.87 

7.17 

5.61 

Nov. 1960 cluster 

8. 57 

4.49 

2. 82 

1.44 

.92 

.66 

. 52 

July 1961 cluster 
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o 

cated at a distance of 1 astronomical unit (1. 496x10 km) from the Sun. Table n(a) pre- 
sents the rad (cj/kg) doses behind various thicknesses of water shield, table n(b) the 
rem doses behind the water shield, table n(c) the rad (cj/kg) doses behind various 
thicknesses of aluminum shield, and table n(d) the rem doses behind the aluminum 
shield. The data presented in this table are input to the computer code MCFLARE. 


Trajectory 

The trajectory selected for the trip is one that requires 556 days for the complete 
mission (and happens to be one that results in near minimum vehicle weight for a depar- 
ture date in 1983). The outward journey from Earth to Mars requires 280 days, then 
there is a 40-day stay at Mars, then a 236-day return trip during which the vehicle ap- 

O 

proaches within 0. 5 astronomical unit (0. 75x10 km) of the Sun. Table HI lists the 


TABLE m. - TRANSFER ELLIPSE PARAMETERS 


Parameter 

Earth to 

Mars to 


Mars 

Earth 

Initial true anomaly, deg 

-24.48 

179. 33 

Final true anomaly, deg 

192. 57 

486. 33 

Eccentricity 

. 1803 

.4687 

Semilatus rectum, km 

1.712X10 8 

1. 099X10 8 


parameters which specify the outbound and return transfer ellipses. These are input to 
the computer code. Trajectory details are discussed in appendix A. 


Discussion of Results 

Figure 1 shows the distribution of dose obtained for a large number of identical 
trips. The figure presents the rem dose distribution behind 20 grams per square centi- 
meter of water shield obtained from a computer run of 40 000 trips. The effect of dis- 
tance on dose was assumed to vary as (1/r) . 

In figure 1(a), A(D), the number of trips on which a dose between D - 1 and D has 
been encountered, is plotted against D. A smoothed curve of this distribution A(D) 
(obtained by averaging seven values centered about each A(D)) is also shown as a dashed 
line in this figure; that is, 

A(D) = A(P - 3) + A (D - 2) + . . . + A(D) + . . . + A(D + 3) 

7 


8 






LOO 



rs 

o 


.98 


.96 


.94 


.92 


.90 


. 88 [ 


/ 




x. 


A 


V 


/ 


X 


V 


A 


V 


. 861 - 


30 


40 50 60 70 80 90 

Dose, D, rem 


100 110 


120 


(b) Cumulative probability of trip encountering dose less than D. 
Figure 1. - Concluded. 


i 

130 


In figure 1(b), p(D) the probability of not exceeding any dose D is shown plotted 
against D 


D 


p(D) = - p A(i) 
M f-t 


i=l 


where M is the number of trips considered in a computer run. For the case of fig- 
ure 1(b), 


P (D) = 


40 000 £ =1 


D 

13 A « 


If Dp is defined as that dose which has a probability p of not being exceeded, then 
from figure 1(b), Dq gg for the group of 40 000 trips, occurs in the dose interval be- 
tween 93 and 94 rem and is selected as the upper bound of the interval, namely, 94 rem. 

The dose probability density function f(D) is approximately equal to A(D)/M. At a 
D q gg of 94 rem, the value of f(D), evaluated by using the smoothed value of A(D), is 
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f(D n QQ ) = 26 - = 0. 00065/rem 

u ° yy 40 000 


and the standard deviation from equation (B6) is a n = 

u 0. 99 

There is a 97. 7 percent confidence that the true dose 


D 0. 99 + 2a D Q 9g ; that is ’ D 0. 99 


is less than 95. 5 rem, 


0. 77 rem. 

Dq 99 is less than 


Table IV presents values of Dq 99 determined by using the computer code for rem 
doses and rad (cj/kg) doses behind various thicknesses (expressed as g/ cm 2 ) of alumi- 
num shield and water shield. These values are presented for both the case where the 

n 

dose received from a flare event varies as (l/r) and for the case where there is no ef- 
fect of distance (i. e. , distance effect, g(r) = 1). Also shown in the table are the values 

of ov. associated with each D n QQ and the corresponding values of 
u 0 o 99 u - yy 

D 0-89 +2<t D 0 .99' 


TABLE IV. - VALUES OF D Q 9g AND CORRESPONDING STANDARD DEVIATION BEHIND VARIOUS THICKNESSES 
OF WATER AND ALUMINUM DETERMINED FROM GROUP OF 40 000 TRIPS 


(a) Rad (cj/kg) dose behind water shield 


(b) Rem dose behind water shield 


Distance 

effect, 

g(r) 


l/r 2 


1 


Quantity 


D 0, 

D 0. 99 + ‘ u D n 


u 0. 99 


0. 99 
» + 2(7 ~ 



Distance | 
effect, 
g(r) 


l/r 2 


Quantity 


Shield thickness, 
g/cm 2 



10 

15 

20 

30 

40 




Dose, 






rem 



D 0. 99 

254 

139 

94 

55 

37 

U 0. 99 

2.0 

1.3 

.77 

.44 

.31 

D 0 99 + 2tJ n 
* U 0 . 99 

258 

141.6 

95. 5 

55.9 

37.6 

D 0. 99 

165 

90 

63 

35 

23 

U 0. 99 

1.7 

1.0 

. 57 

i .25 

.2 

D 0 99 + 2(7 d 
* U 0. 99 

168.4 

92 

64. 1 

35. 5 

23.4 


(c) Rad (cj/kg) dose behind aluminum shield 


(d) Rem dose behind aluminum shield 


Quantity 



Shield thickness, 






g/cm 2 





10 

15 

20 

30 

40 

50 

60 





Dose, 







rad (or cj/kg) 



D 0. 99 

327 

178 

116 

63 

41 

29 

21 

X.99 

2.9 

1. 3 

1.2 

. 54 

.33 

.25 

. 18 

D 0-99 + 2a D 0 99 

332.8 

180. 6 

118.4 

64. 1 

41.7 

29.5 

21.4 

D 0. 99 

217 

115 

75 

41 

26 

18 

13 

X.99 

2.0 

1.2 

.67 

.29 

.2 

. 14 

. 10 

D 0-99 +2 X.99 

221 

117.4 

76.3 

41.6 

26.4 

18.3 

13.2 


Distance 

effect, 

B(r) 


l/r 2 


1 


Quantity 



Shield thickness, 
g/cm 2 




10 

15 

20 

30 

40 

50 

60 





Dose, 








rem 




D 0. 99 

376 

213 

146 

87 

62 

48 

39 

^0. 99 

3.3 

1.8 

1.4 

.77 

. 54 

.42 

.31 

D 0 nq + 2a^ 
u. aa d q> 9Q 

382.6 

216.3 

148. 8 

88. 5 

63. 1 

48.8 

39.6 

D 0. 99 

250 

138 

94 

58 

40 

31 

25 

Qd 0.99 

2.2 

1. 5 

.91 

. 59 

.31 

.23 

.22 

D 0-99 +2a D 0 .99 

254.4 

141 

95.8 

59.2 

40.6 

31. 5 

25.4 
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Aluminum shield thi( 


(c) Rad IcJ/kg) dose behind aluminum shield. 


Figure 2. - 



(d) Rem dose behind aluminum shield. 




In figure 2 the values of Dq gg + 2 ct d are plotted against shield thickness. In 

0. 99 

figure 2(a), the rad (cj/kg) dose is plotted against water shield thickness for cases of 

p 

g(r) = 1/r and g(r) = 1. Figure 2(b) is a similar plot of the rem dose behind the water 
shield, figure 2(c) is a similar plot of rad (cj/kg) behind the aluminum shield, and fig- 
ure 2(d) a similar plot of rem dose behind the aluminum shield. 

These curves indicate how much more effective water is as a shield material than 
aluminum for shielding against solar flare protons on the basis of grams per square 
centimeter of material necessary to maintain a given dose level inside the crew com- 

p 

partment. Also shown for this mission is that a 1/r distance effect on dose received 
from a flare can have an appreciable effect on shield requirements when the mission 

Q 

trajectory brings the vehicle to within 0. 5 astronomical unit (0. 75x10 km) of the Sun. 
This result indicates a need for information regarding the effect of position from the Sun 
on radiation encountered from a flare. 

COMPUTER PROGRAM 

Complete data input instructions for the computer program MCFLARE are pre- 
sented in this section. In its present form, MCFLARE will accept as many as 25 differ- 
ent flare types with up to seven different shield configurations. Doses are tallied in as 
many as 900 unit dose bins. With these limits, the program does not occupy the maxi- 
mum 32 K storage available in the core. The user can increase the limits of any or all 
these arrays. The version of the code presented herein is operational on the Lewis 
Research Center IBM 7094-11/7044 computer system. Execution times for the examples 
discussed previously in this report were about 2 minutes per group of 40 000 trips. 
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Data Input Flow Chart 



Card column Variable 

1-5 MTHOUS 

6-10 NTHOUS 

11-15 IXP 

16-20 NTYPE 

V 

Go to next page 


Read XYZ. 

XYZ is the last random number used in a 
previous run. This number is punched 
out at the end of each case. If calculation 
is being started without a calculated XYZ, 
XYZ = 000000000001 should be used. 

XYZ is a 12 -place octal integer, and 
punched in card columns 1 to 12. 

Read 3 cards. Format (3(12A6)) 

These cards contain any alphanumeric in- 
formation in columns 1 to 72. This infor- 
mation is used for titling output. 

Read control card. Format (615, 2E10. 4) 


MTHOUS is the number of thousands of trips 
considered. An intermediate printout is 
made after each NTHOUS thousand. If 
NTHOUS = 0, no intermediate printout 
will be made. 

IXP 

Exponent in (1/r) dependence of dose 
with distance from Sun. If IXP = 0, there 
is no effect of distance on dose. (IXP is 
called a in this report. ) 

Number of different particular flare events 
considered (NTYPE < 2 5) 


* 



Symbol means read a card or set of cards as described. 
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Card column Variable 


21-25 NSHELD 


26-30 NBIN 


31-40 TREF 


41-50 TOTALT 



Variable 

UNITS 


SHIELD j 


shield 2 


Number of different shields considered 
(NSHELD s 7) 

Number of dose bins of unit width to be used 
(NBIN < 900); (A final bin extending to in- 
finity is used by the program to maintain 
a correct tally in each bin. ) 

The duration of the solar active period is 
given in days. Each flare event is as- 
sumed to occur with expected frequency of 
once every TREF days. 

The length of the mission is given in days. 
This parameter is used only if IXP = 0 
because no radius (trajectory) is calcu- 
lated in this case. If IXP * 0, TOTALT 
is calculated from input trajectory param- 
eters. If IXP = 0, TOTALT must be 
given here. 

Read output titling card. Format (8A6) 


The dose units used in this problem, for 
example, RADS 

A name for the first shield, for example, 
10GMAL for 10 grams per square centi- 
meter of aluminum 

A name for the second shield 

(NSHELD values required) 
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c 


5 


Card column 


1-9 


10-18 

19-27 


Read dose table. Format (7E9. 4) 

For each particular flare, read a card list- 
ing the dose at 1 astronomical unit 
(1. 496X10 8 km) with each NSHELD shield 
in place. One card is required for each 
particular flare used. NSHELD values 
are required per card. 


Dose with shield 1 for this particular flare 
Dose with shield 2 for this particular flare 
Dose with shield 3 for this particular flare 



TDUR(l) 

TDUR(2) 


Read flare duration card(s) Format (8E9. 4) 
For each particular flare event read in the 
flare duration in days. If clustered flares 
are considered to be single events, an ap- 
propriate time must be allotted. NTYPE 
values are required. 


Time duration of particular flare 1 
Time duration of particular flare 2 
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15-22 THETAFj 

23-30 E* 


Punch out last 
random number 

40-47 THETA 2 


48-55 THETAF 2 


56-63 E 2 


Go back to read 
card 2 and con- 
tinue with next 
case 


31-39 P x 



Read mission trajectory card. 
Format (F6. 0, 3F8. 0, E9. 4, 

3F8. 0, E9. 4) 

Number of days stay time at planet 

True anomaly (angle from perihe- 
lion) at the start of first transfer 
ellipse is in degrees. 

True anomaly at the end of first 
transfer ellipse is in degrees. 

Eccentricity of first transfer ellipse 

Semilatus rectum of first transfer 
ellipse is given in kilometers. 

True anomaly at the start of second 
transfer ellipse is given in 
degrees. 

True anomaly at the end of second 
transfer ellipse is given in 
degrees. 

Eccentricity of second transfer 
ellipse 


a 


64-72 P 2 Semilatus rectum of second trans- 

fer ellipse is given in kilometers. 


(NOTE: THETAF must be greater 
than THETA. ) 
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Sample Problem Input 

To illustrate data input to the computer code, consider the case with water shield- 

p 

ing, rem dose, 1/r distance dependence of dose, and the 556-day mission described in 
the section ILLUSTRATIVE EXAMPLE. The complete set of data input is illustrated in 
table V. 

Parameters on the control card include 

MTHOUS = 40 40 000 Trips will be considered 

NTHOUS =0 No intermediate printout will be made 

IXP =2 1/r^ Spatial dependence 

NTYPE = 13 13 Flare events will be considered, as indicated in table n 

NSHELD =5 5 Shields will be used 

NBIN = 500 As many as 500 dose bins will be used 

TREF = 2000. 0 The observed flares were assumed representative of a 

2000. 0-day solar active period 

The titling card includes 

UNITS = REMS Dose tables input in rems 

SHIELD i = 10 H20 10 g/cm^ water is first shield 

SHIELDg = 15 H20 15 g/cm^ water is second shield 

SHIELD^ = 20 H20 20 g/cm^ water is third shield 

SHIELD^ = 30 H20 30 g/cm^ water is fourth shield 

SHIELDg = 40 H20 40 g/cm^ water is fifth shield 

The dose table is taken directly from table 11(b). 

The duration time of each event, listed on the next cards, was taken to be 2. 0 days 
for all single events, 8. 0 days for the clustered events of July 1959 and July 1961, 10. 0 
days for the events of November 1960, and 12 days for the August 1958 events. The tra- 
jectory information on the next card is taken from table m. The hold time at the planet, 
THOLD, was taken to be 40. 0 days. Note the subroutine RADIUS requires that the final 
true anomaly be greater than the initial true anomaly; thus, in the second transfer el- 
lipse, the final true anomaly is given as 486. 33 rather than 126. 33 degrees. In the first 
transfer ellipse, the true anomaly ranges from -24. 48 to 192. 57 degrees. This transfer 
ellipse could also have been specified as 0 = 6 + 360° starting at 335. 52 and ending at 
552. 57 degrees. 
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TABLE V. - DATA INPUT FOR SAMPLE CASE 


TITLE 



MCFLAHE 




PROJECT NUMBER 

ANALYST 

SHEET OF 

STATEMENT 

NUMBER 

i 

FORTRAN STATEMENT 

IDENTIFICATION 

1 2 3 4 5 

D 

7 8 9 10 11 12 

IS 14 IS 16 17 IB 

19 20 

21 22 23 24 

25 26 27 28 29 30 

31 32 33 34 35 36 

37 38 39 40 

41 42 

43 44 45 46 47 48 

49 50 51 52 53 54 

55 56 57 58 59 60 

61 62 63 64 65 66 

67 68 69 70 71 72 

73 74 75 76 77 78 79 80 

0,0, 0,0,0 

0 

0000.0.1 
















,P,R,Op,L 

EM. . . . 

. . .5.5 .6 

-DAY. M 

A,R£, 

MI 

SSION. 

, , P,.P 

A|D, ,A,P 

PR,OA,0|fi 


Wf#) 



WATER, 

SHIELD 

S 


, 4^6- DAY 

a^s 

STOPOVER 





C,0f4,M^EN T,2. 



REM. DOSES. . . 



1 / R*.*.g. 

SPAT 

i 

I.AL DEPENDENCE 






f! l o,M l M,F..WT'3 

f 4,0 


,0, 4-. 

, P, , , 
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_ -2 
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1.09 


0 
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Sample Problem Output 


The computer output generated from execution of the sample problem detailed in the 
previous section is shown in table VI. The first page of computer output is a printout of 
the input information of cards 2 to 6. If radius as a function of trip time is calculated 
(IXP 4 - 0), the next page of computer output is the trajectory information, namely, a 
listing of mission parameters and a 100-point table listing time elapsed in the trip and 

O 

the corresponding distance to the Sun in astronomical units (1. 496x10 km) at that time. 

The next page of computer output is a table of the number of trips observed in this 
group having encountered the particular flares n times (n = 0, 1, 2, . . . ). For ex- 
ample, the first input flare event did not occur on 30 586 of the executed 40 000 trips, 
occurred once on 8223 trips, occurred twice on 1092 trips, occurred three times on 92 
trips, and occurred four times on 7 trips; all 40 000 trips were accounted for. The last 
column of this table indicates the number of trips observed which encountered n flares, 
independent of which particular flare occurred. For example, in table VI, 1074 trips 
encountered no flares at all, while 2 trips experienced a total of 13 flares. For com- 
parison, the expected number per trip of each particular flare type (equal to the ratio of 
total trip time to reference time) and of any flare type (equal to the product of this ratio 
and the number of flare types considered) are calculated. Then, for an assumed 
Poisson distribution, the probability per trip of n flares (n = 0, 1, 2, . . . , 10) occurring 
are calculated and printed out. It should be pointed out here that the observed frequency 
tables should agree with this calculated table only in the case where the time allowed for 
a flare event to occur goes to zero. In this case, the problem reduces to a constant- 
rate process, one that obeys a Poisson distribution. 

Finally, the tabulation of number of trips observed to have a dose lying in a unit 
dose bin against dose is printed out for each shield configuration. The first column is 
the upper limit of the dose bin of unit width; in the next column, for each shield is the 
number of trips which encountered a dose lying in this unit dose bin followed by the per- 
centage of trips observed to have a dose less than or equal to this upper limit dose. 
From this list (or a smoothed one) the dose which will not be exceeded p percent of the 
time, D may be ascertained for each shield. 

Jr 
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TABLE VI. - COMPUTER OUTPUT OF SAMPLE PROBLEM 


SAMPLE PRCBLEM 
WATER SHIELDS 
REM DOSES 


556 DAY MARS MISSIS 0.5 AU APPROACH 

AC DAY MARS STOPOVER 
1/R **2 SPATIAL DEPENDENCE 


THIS CALCULATION CONSIDERS 13 FLARES, AND 5 SHIELDS, NAMELY 10 H20 15 H20 20 H20 30 H20 AO 

DC SE CALCULATICN USES UP TO 500 UN IT BINS 
AC, 000 TRIPS TO BE RLN WITH PPINTCUT AFTER EACH AO, 000 
DOSES EASED ON l./R** 2 


FLARE INPUT DATA 


DOSES IN *REMS* 


FLARE OCCURRENCES OF DURATION COSE FROM ITH TYPE FLARE 


TYPE. 

ITH TYPE FLARE 

OF ITH TYPE 

AT 1 

AU 

I 

IN 2CCC .0 DAYS 

FLARE, DAYS 

10 

H2C 

1 

1 

C.20CCE 

01 

0.A210E 

C2 

2 

1 

C.20CCE 

01 

0.1C7CE 

01 

3 

1 

C.20CCE 

01 

C.A300E 

CO 

A 

1 

C.20COE 

01 

C.1200E 

Cl 

5 

1 

C.2CCCE 

Cl 

0.U20E 

01 

6 

1 

C.20CCE 

Cl 

0. 12606 

01 

7 

1 

C.I200E 

02 

C.1030E 

01 

8 

1 

C.20CCE 

Cl 

0. U30E 

C2 

9 

1 

C.20CCF 

Cl 

C.6Q00E 

CO 

10 

1 

C.800CE 

01 

0.A959E 

C2 

11 

1 

C.20C0E 

01 

C.8AOOE 

00 

12 

1 

C. 10CCE 

C2 

C.A93AE 

02 

13 

1 

C.80CCE 

01 

0.5A80E 

Cl 


WITH JTH SHIELD- 


15 

H20 

20 

H20 

30 H20 

AC H20 

0.2776E 

02 

0.2079E 

02 

0.1303E 02 

C.903CE 01 

C.380CE 

CO 

0. 190CE 

OC 

C.700CE-C1 

C.AOOCE-Ol 

0.150C6 

00 

0 • 70C0E- 

-01 

0.30COE-C1 

0. 10CCE-01 

C.650CE 

00 

0.A2CCE 

00 

0.2100E CC 

C.120CE 00 

C.AOOCE 

00 

0. 190 C E 

00 

C.7000E-0 1 

C.A0CCE-C1 

0.A50CE 

00 

0.22CCE 

OC 

0.8000E-0 l 

C.AOCCE-Ol 

0.350CE 

00 

0. 17CCE 

00 

C.7000E-01 

O.AOOOE-Ol 

0.A930E 

01 

0.2670E 

01 

0.1090E 01 

C.570CE 00 

0.2000E 

00 

0. 100CE 

00 

O.AOCOE-Cl 

C.20CCE-C1 

0.2368E 

02 

0.1376E 

02 

0.611 CE 01 

0.336CE 01 

O.A600E 

00 

0.29C0E 

00 

0.1500E CO 

C.90C0E-01 

0.26A9E 

02 

0. 1679E 

02 

0.8300E Cl 

0.A85CE Cl 

0.2670E 

Cl 

0. 1570E 

01 

C.7100E 00 

C.390CE 00 


MISSION PARAMETERS 

TOTAL TRIP TIME. DAYS 556.523C 


EARTH TO PLANET, 
SEMI-LATUS RECTLM C.1712CE C9 
ECCENTRICITY C.10C3O 

TRUE ANOMALY AT START -2A.AeO 
TRUE ANOMALY AT END 192. 57C 
PERIOD. CAYS A69.S112 

SEMI-MAJOR AXIS 0.17695E C9 


TC EARTH 
C.IC99CE C9 
C.A687C 
179. 33C 
A86.330 
333.6709 
C.1AC8AE 09 


KM. 

DEGREES 

DEGREES 

KM 


ELAPSED TIME 

(DAYS) 

vs m 

iSTANCE FROM SUN 

{ AU) 








1 

0 . 

C.983 

21 

82.915 

1.C61 

A 1 

207.336 

1.363 

61 

429.898 

0.942 

81 

490.659 

2 

A • 0 A 2 

C.S79 

22 

fi 7. 7C0 

1 .074 

42 

215. 195 

1.372 

62 

436.113 

0.889 

82 

492.535 

3 

8.05C 

C.S75 

23 

92.6C3 

1.087 

A3 

223. 156 

1.381 

63 

A A 1 • 652 

C.840 

83 

494.442 

A 

12.033 

C.972 

2A 

97.631 

1.101 

AA 

231.203 

1.387 

64 

446.600 

0.794 

84 

496.396 

5 

15.996 

C.971 

25 

1 02. 793 

1.116 

45 

239.316 

1.392 

65 

A51.C36 

0.753 

05 

490.412 

6 

19.952 

C.S70 

26 

108. C95 

1.L31 

46 

247.475 

1.395 

66 

A55.03C 

0.715 

86 

5CC.510 

7 

23.902 

C.S70 

27 

113. 5A5 

1.147 

A7 

255.658 

1.396 

67 

458. 6A5 

0.682 

87 

5C2.7C8 

8 

27.855 

C.971 

28 

119.150 

1.163 

48 

263.845 

1.396 

68 

461.937 

0.651 

88 

5C5.C29 

9 

31.819 

C.972 

29 

12 A. 9 1 5 

1.180 

49 

272.014 

1.393 

69 

464.953 

0.625 

89 

507.498 

10 

35.801 

C.975 

3 C 

13C.8A5 

1.197 

50 

280. 1AA 

1.389 

70 

467.735 

0.601 

90 

5LC.143 

11 

39.8C8 

C.978 

31 

136. 9A6 

1 .2 1 A 

51 

320. 1AA 

1.383 

71 

470.319 

0.58C 

91 

512.997 

12 

A3 • 8 A 6 

C.983 

32 

143.220 

1.231 

52 

33A.289 

1 .377 

72 

472.736 

0.562 

92 

516.098 

13 

A 7.927 

C.988 

33 

1A9.669 

1.2A8 

53 

3A8. 173 

1.357 

73 

475.014 

0.547 

93 

519.490 

1A 

52 - 0 5 A 

C .99 A 

3 A 

156.295 

1.26A 

54 

361.529 

1.325 

74 

477.177 

0.534 

94 

523.223 

15 

56.236 

1.CC1 

35 

163. C97 

1.281 

55 

374.143 

1.282 

75 

479.247 

0.523 

95 

527.354 

16 

60 .A 8 1 

1.009 

36 

17C.073 

1.296 

56 

385.862 

1.231 

76 

48 1 .241 

0.514 

96 

531.948 

17 

6A.795 

1.C18 

37 

177.217 

1.312 

57 

396.602 

1.175 

77 

483.178 

0.508 

97 

537. 081 

18 

69-188 

1.C27 

38 

1 8 A. 52 5 

1.326 

58 

A06.338 

1.116 

78 

485.073 

0.503 

98 

542.831 

19 

73.667 

.1 • C 3 R 

39 

191.980 

1.339 

59 

415. C92 

1.057 

79 

486.942 

0.501 

99 

549.286 

20 

78.2AC 

1 • C49 

AO 

159.595 

1.352 

60 

422.920 

0.998 

80 

408.799 

C.50C 

ICC 

556.533 


H20 


C.502 
0.505 
0.510 
0.517 
0.526 
C.538 
0.552 
0.568 
C. 587 
0.6C8 
0.633 
0.661 
0.692 
C .727 
C. 766 
C.809 
0.855 
C .906 
C.960 
1.017 
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TABLE VI. - Continued. COMPUTER OUTPUT OF SAMPLE PROBLEM 


THE NUMBER OF TRIPS IN WHICH I FLARES CF TVPE J CCCUR — 

I THE LAST COLUMN IS THE NUMBER OF TRIPS IN WHICH I FLARES OF ANY TYPE OCCUR) 


I 

J = 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

0 

30586 

3C433 

30585 

30556 

30554 

3C445 

30573 

30501 

3 04 8 A 

30446 

30459 

30540 

30467 

1074 

1 

8223 

63 16 

8155 

8273 

8183 

8341 

8277 

8182 

8302 

8316 

8390 

8252 

8338 

4016 

2 

1092 

1122 

1175 

1096 

1164 

1097 

1057 

1193 

1113 

1129 

1075 

1102 

1092 

7366 

3 

92 

12C 

77 

7 C 

91 

113 

89 

1 19 

87 

105 

72 

97 

98 

8703 

A 

7 

e 

e 

5 

8 

4 

4 

A 

1 1 

4 

4 

9 

5 

7902 

5 

0 

l 

0 

C 

0 

0 

0 

1 

3 

0 

0 

0 

0 

5412 

6 

0 

c 

c 

C 

0 

0 

c 

0 

c 

0 

C 

0 

C 

3C89 

7 

0 

0 

c 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1454 

8 

0 

c 

c 

c 

0 

0 

0 

c 

0 

0 

0 

0 

0 

643 

9 

0 

c 

0 

c 

0 

0 

0 

c 

0 

0 

0 

0 

c 

233 

IC 

c 

c 

c 

c 

0 

0 

0 

0 

C 

0 

C 

0 

c 

75 

11 

c 

c 

c 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

25 

12 

c 

c 

c 

c 

0 

0 

0 

c 

C 

0 

0 

0 

0 

6 

13 

c 

c 

c 

c 

0 

0 

0 

0 

C 

0 

0 

c 

0 

2 

SUMS 

4C0GQ 

4CCCC 

4CCCC 

40000 

40000 

40000 

40000 

40000 

4000 C 

40000 

4 00 0 C 

AO 000 

40000 

40000 


THE EXPECTED NO. OF EACH TYPE IS ... 0.278 

THE EXPECTED NO. OF ANY TYPE IS ... 3.617 


CALCULATED PROBABLE OCCURRENCES PER TRIP USING PCISSON DISTRIBUTION 


NUMBER 

EACH TYPE 

ANY TYPE 

0 

7.57C55E-01 

2.665C7E-C2 

l 

2. 1C67AE— Cl 

9.7131 3E— C 2 

2 

2 .531 18E— 02 

1.75685E-C1 

3 

2.718836-03 

2.11844E-C1 

4 

1.8914CE-C4 

1.915 8 5E— C 1 

5 

1.C5263E-C5 

1.3661CE-CI 

6 

4.88 1 84 E— 0 7 

8.35656E-C2 

7 

1 . 9 406 5 E— 0 8 

4.31871 E— C 2 

8 

6.75C21E— 10 

1.95285E-C2 

9 

2.CE7C6E-11 

7.8AS29E-C3 

10 

5.8076CE-13 

2.835A5E-C3 
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TABLE VI. - Continued. COMPUTER OUTPUT OF SAMPLE PROBLEM 


COSE OISTRIBLTICN 

SHIELD— 1C H20 SHIELD-15 H20 SHIELD-20 H20 SHI ELD- 3C H20 SHIELD-40 H20 SHIELD- SH1EL0- 



NO. 

PERCENT 

NO. 

PERCENT 

NO. 

PERCENT 

NO. 

PERCENT 

NO. 

PERCENT 

NO. PERCENT 

NO. PERCENT 

OOSE *R EM S* 

IN 

LESSTHAN 

IN 

LESSTHAN 

IN 

LESSTHAN 

IN 

LESSTHAN 

IN I 

LESSTHAN 

IN LESSTHAN 

IN LESSTHAN 

UPPER L IM IT 

BIN 

LIMIT 

BIN 

LIMIT 

eiN 

LIMIT 

BIN 

LIMIT 

BIN 

LIMIT 

BIN LIMIT 

BIN LIMIT 

1 

3147 

7.667 

6764 

16.910 

9296 

23.240 

13489 

33.722 

160C6 

40.C15 



2 

2592 

14.347 

3C35 

24.497 

3485 

31.952 

2701 

40.475 

2110 

45.29C 



3 

168 7 

19.C65 

1585 

29.460 

1906 

36.717 

795 

42.462 

3237 

53.382 



4 

1339 

22.412 

1565 

33.372 

1C9 5 

39.465 

1675 

46.650 

2400 

59.382 



5 

11 38 

25.257 

549 

35.745 

555 

40.852 

2143 

52.0C7 

2101 

64.635 



6 

963 

27.665 

815 

37.782 

380 

41.802 

1566 

55.922 

2311 

70.412 



7 

879 

29.862 

586 

39.247 

280 

42.502 

1708 

60.192 

1459 

74.060 



8 

764 

31.772 

356 

40.137 

856 

44.642 

1700 

64.442 

1249 

77.182 



9 

572 

33.2C2 

3C5 

4C.900 

967 

47.060 

1407 

67.960 

1136 

80.C22 



1C 

500 

34.452 

219 

41.447 

1 19 C 

5C.035 

1141 

70.812 

1255 

83. 160 



11 

396 

35.442 

214 

41.982 

1185 

52.997 

951 

73.19C 

852 

85.29C 



12 

391 

36.42C 

146 

42.347 

1377 

56.440 

870 

75.385 

659 

86.937 



13 

346 

37.265 

458 

43.492 

1023 

58.997 

779 

77.332 

604 

88.447 



14 

316 

38.C75 

618 

45.C37 

851 

61.125 

923 

79.64C 

546 

89.812 



15 

276 

38.765 

1 C52 

47.767 

843 

63.232 

769 

81.562 

407 

9C.830 



16 

221 

39.317 

1112 

50.547 

737 

65.075 

629 

03.135 

394 

91.815 



17 

184 

39.777 

911 

52. 825 

642 

66.680 

497 

84.377 

331 

92.642 



18 

167 

40.195 

758 

54.820 

660 

60.33C 

480 

85.577 

310 

93.417 



19 

152 

40.575 

753 

56.7C2 

705 

7C.092 

425 

86.64C 

267 

94.C85 



20 

108 

40.645 

688 

58.422 

603 

71.600 

362 

87.545 

262 

94.740 



21 

122 

41. 15C 

59 C 

59.897 

548 

72.970 

336 

88.385 

185 

95.202 



22 

182 

41.6C5 

5C7 

61.165 

570 

74.395 

320 

89.185 

190 

95.677 



23 

291 

42- 332 

438 

62.260 

628 

75.965 

292 

89.915 

147 

96.045 



24 

276 

43. C22 

433 

63.342 

506 

77.23C 

296 

90.655 

138 

96.390 



25 

273 

43. 7C5 

440 

64.442 

443 

78.337 

282 

91.360 

119 

96.687 



26 

399 

44.7C2 

502 

65.697 

391 

79.315 

219 

91.907 

91 

96.915 



27 

541 

46.C55 

4 73 

66.880 

374 

80.250 

215 

92.445 

97 

97. 157 



28 

572 

47.485 

406 

68.095 

328 

8 1.070 

184 

92.905 

82 

97.362 



29 

566 

48.9CC 

559 

69.492 

332 

81.90C 

188 

93. 375 

92 

97.592 



30 

515 

5C. 167 

554 

7C.877 

308 

82.670 

165 

93.787 

65 

97.755 



31 

454 

51.322 

5C8 

72.147 

261 

83.322 

174 

94.222 

80 

97.955 



32 

454 

52.457 

368 

73.117 

305 

84.085 

147 

94.590 

75 

98. 142 



33 

43 1 

53.535 

363 

74.C25 

271 

84.762 

164 

95.0CC 

59 

98.29C 



34 

399 

54.532 

348 

74.895 

271 

85.44C 

168 

95.420 

64 

98.450 



35 

402 

55.537 

316 

75.685 

219 

85.987 

107 

95.687 

84 

98.66C 



36 

332 

56.367 

281 

76.387 

199 

86.485 

103 

95.945 

82 

98.065 



37 

357 

57.26C 

268 

77.C57 

192 

86.965 

93 

96.177 

84 

99.C75 



38 

344 

58.120 

2 €5 

77.770 

200 

87.465 

111 

96.455 

44 

99.185 



39 

312 

58.9CC 

255 

78.407 

205 

87.977 

87 

96.672 

59 

99.332 



40 

329 

55.722 

268 

79 .077 

178 

88.422 

74 

96.857 

39 

99.430 



41 

293 

60.455 

232 

79.657 

140 

88.772 

71 

97.035 

38 

99.525 



42 

253 

61.C87 

269 

8C.330 

174 

89.207 

79 

97.232 

23 

99.582 



43 

262 

61.742 

256 

80.970 

149 

89.580 

62 

97.387 

25 

99.645 



44 

262 

62.397 

217 

81.512 

144 

89.940 

55 

97.525 

11 

99.672 



45 

332 

63. 227 

231 

82.C90 

134 

90.275 

49 

97.647 

20 

99.722 



46 

3C9 

64.CCC 

206 

82.6C5 

135 

9C.612 

73 

97 . 83 C 

20 

99.772 



47 

269 

64.672 

217 

03.147 

115 

9C.90C 

58 

97.975 

15 

99.810 



48 

269 

65.345 

182 

83.602 

12 1 

9 1.202 

53 

98. 107 

7 

99.827 



49 

241 

65.547 

162 

84.0C7 

130 

91.527 

45 

58.22C 

11 

99.855 



50 

262 

66.602 

173 

84.440 

125 

91.840 

50 

98.345 

12 

99.805 



51 

244 

67.212 

165 

84. 862 

104 

92.100 

66 

98.51C 

5 

99.897 



52 

299 

67.960 

169 

85.285 

111 

92.377 

70 

98.685 

6 

99.912 



53 

338 

60.8C5 

135 

85.622 

110 

92.652 

63 

98.842 

3 

99.920 
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TABLE VI. - Continued. COMPUTER OUTPUT OF SAMPLE PROBLEM 


54 

29 8 

6S.55C 

136 

85.962 

122 

92.960 

47 

98.960 

5 

99.932 

55 

318 

7C.345 

14C 

86.312 

111 

93.237 

36 

99.050 

3 

99.940 

56 

256 

70*585 

148 

€6.682 

120 

93.537 

36 

99.140 

4 

99.950 

57 

248 

71*605 

11C 

86.557 

103 

93.795 

37 

99.232 

l 

99.952 

58 

236 

72.195 

140 

87.307 

89 

94.017 

35 

99.320 

3 

99.960 

59 

2C9 

72.71? 

143 

87.665 

77 

94.210 

34 

99.405 

l 

99.962 

60 

215 

73.255 

110 

87.540 

96 

94.45C 

24 

99.465 

3 

99.970 

61 

161 

73.(57 

113 

88.222 

83 

94.657 

14 

99.500 

l 

99.972 

62 

186 

74.122 

116 

88.512 

65 

94.82C 

16 

99.540 

1 

99.975 

63 

193 

74.6C5 

101 

88.765 

75 

95.007 

15 

99.577 

0 

99.975 

64 

165 

75.C17 

95 

89.C02 

70 

95.182 

17 

99.620 

1 

99.977 

65 

164 

75.427 

• ICO 

89.252 

09 

95.405 

19 

99.667 

2 

99.982 

66 

162 

75. €32 

95 

89.500 

84 

95.615 

6 

99.682 

2 

99.987 

67 

137 

76.175 

59 

89.747 

84 

95.825 

14 

99.717 

0 

99.987 

68 

152 

76.555 

89 

89.970 

68 

95.995 

13 

99.75C 

0 

99.587 

69 

160 

76.555 

93 

90.202 

77 

96. 187 

6 

99.765 

0 

99.987 

70 

155 

77.342 

89 

90.425 

70 

96.362 

9 

99.787 

1 

99.990 

71 

162 

77.747 

86 

90-640 

63 

96.520 

6 

99.802 

0 

99.990 

72 

151 

78.125 

86 

50.855 

52 

96.650 

6 

99.817 

1 

99.992 

73 

151 

78.5C2 

89 

51.077 

62 

96.805 

8 

99.837 

L 

99.995 

74 

147 

78.67C 

54 

91.312 

47 

96.922 

5 

99.850 

0 

99.995 

75 

157 

79.262 

58 

91.457 

48 

97.042 

4 

99.860 

0 

99.995 

76 

124 

79.572 

76 

91.647 

58 

97. 187 

8 

99.880 

1 

99.997 

77 

151 

75.55C 

71 

91.825 

59 

97.335 

3 

99.887 

0 

99.997 

78 

133 

SC. 282 

72 

92 • C 05 

48 

97.455 

6 

99.902 

0 

99.997 

79 

161 

8C.685 

61 

92. 157 

44 

97.565 

3 

99.910 

1 

ICC. 000 

80 

141 

81. €37 

84 

92.367 

51 

97.692 

3 

99.917 

G 

ICO. 000 

81 

148 

01.4C7 

73 

52.550 

43 

97.80C 

2 

99.922 

C 

ICO. 000 

82 

121 

61.71C 

71 

92.727 

41 

97.902 

2 

99.927 

c 

ICO. 000 

83 

118 

82.CC5 

57 

92.870 

67 

98.070 

I 

99.930 

c 

ICO. 000 

84 

125 

82.31? 

62 

53.025 

54 

98.205 

7 

99.947 

0 

100.000 

85 

109 

82.59C 

61 

53.177 

46 

98.320 

4 

99.957 

0 

ICO. 000 

86 

121 

82. 892 

58 

93.322 

33 

98.402 

l 

99.960 

c 

ICO. coo 

87 

1C 7 

83. 16C 

58 

53.467 

30 

98.477 

0 

99.960 

c 

ICO. 000 

88 

104 

83.420 

77 

93.660 

31 

98.555 

0 

99.960 

c 

ICO. 000 

89 

1 1C 

83.655 

(8 

93.830 

37 

98.647 

1 

99.962 

c 

ICO. 000 

90 

96 

83.935 

59 

53.977 

24 

98.707 

1 

99.965 

0 

ICO. 000 

91 

93 

84. 167 

69 

54.150 

27 

98.775 

2 

99.970 

0 

ICO. coo 

92 

86 

64.382 

71 

54.327 

28 

98.845 

2 

99 .975 

0 

ICO. 000 

93 

72 

84*562 

52 

94.457 

35 

98.932 

1 

99.977 

0 

ICO. 000 

94 

94 

84.797 

73 

94 • 64 C 

38 

99.027 

1 

99.980 

c 

100.000 

95 

£ 8 

85. Cl? 

66 

94.005 

18 

99.072 

2 

99.985 

0 

ICO. 000 

96 

85 

85.23C 

76 

54.995 

16 

99.112 

C 

99.985 

c 

100.000 

97 

75 

85.417 

53 

95.127 

15 

99. 15C 

1 

99.987 

c 

100.000 

98 

80 

85.(17 

52 

95.257 

12 

99. 180 

C 

99.987 

c 

1 CO. 000 

99 

73 

as.ecc 

53 

55.390 

21 

99.232 

0 

99.987 

c 

ICO. 000 

100 

ICO 

86.C50 

49 

55.512 

16 

99.272 

0 

99.987 

c 

100.000 

101 

81 

f 6. 252 

39 

95.610 

11 

99.300 

0 

99.587 

c 

ICO. 000 

102 

76 

86.442 

48 

55.73C 

20 

99.350 

0 

99.987 

c 

ICO. 000 

103 

84 

86.(52 

51 

55.857 

17 

99.392 

l 

99.990 

c 

ICO. 000 

104 

72 

€6.832 

50 

95.982 

5 

99.405 

0 

99.990 

c 

100.000 

105 

68 

87.CC2 

53 

96.115 

16 

99.445 

0 

99.990 

c 

ICO. 000 

106 

70 

87. 177 

57 

96.257 

13 

99.477 

0 

99.990 

c 

100 . coo 

107 

68 

67.347 

60 

96.407 

14 

99.512 

1 

99.992 

c 

100.000 

108 

60 

€7.497 

(2 

96.562 

15 

99.550 

0 

99.992 

c 

100.000 

109 

69 

67.67C 

56 

56.702 

8 

99.570 

0 

99.992 

c 

ICC. 000 

1 10 

82 

€7.675 

51 

56.830 

9 

99.592 

l 

99.995 

c 

ICC. 000 

ill 

42 

€7.580 

68 

57. COO 

7 

99.610 

0 

99.995 

0 

100.000 

112 

58 

88.125 

50 

97. 125 

1C 

99.635 

0 

99.995 

0 

ICO. 000 

113 

57 

88.267 

39 

97.222 

14 

99.670 

1 

99.997 

c 

100.000 
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114 

66 

88.432 

38 

97.317 

6 

99.685 

0 

59.997 

c 

IOC. ooo 

115 

66 

88.597 

35 

97.405 

7 

99.702 

0 

99.997 

0 

ICO. ooo 

1 16 

69 

88.770 

33 

97.487 

3 

99. 71 C 

c 

99.997 

c 

ICO. ooo 

117 

60 

08.920 

38 

97.582 

6 

99.725 

0 

59.997 

c 

100.000 

118 

62 

89.C75 

38 

97.677 

3 

99.732 

I 

ICO. 000 

c 

ICO. ooc 

1 19 

55 

89.212 

31 

97.755 

8 

99.752 

c 

ICO. 000 

0 

ICO. ooo 

120 

40 

89.222 

33 

97.837 

4 

99.762 

c 

ICO. 000 

c 

ICO. ooo 

121 

51 

89.460 

52 

97.967 

5 

99.775 

0 

100.000 

0 

ICO. coo 

122 

49 

89.582 

32 

96.C47 

4 

99. 785 

0 

10C. 000 

0 

ICO. ooo 

123 

55 

89.720 

25 

98. U0 

2 

99. 79C 

0 

ICC. 000 

0 

ICO. OOO 

124 

45 

89.622 

3 L 

98. 187 

3 

99.797 

0 

100.000 

c 

ICC. OOO 

125 

55 

89.970 

25 

90.250 

1 

99. 800 

0 

ICC. 000 

c 

100.000 

126 

40 

9C.C90 

27 

98.317 

4 

99.810 

0 

100.000 

0 

ICO. ooo 

127 

56 

9C.23C 

34 

98.402 

3 

99.817 

0 

ICC. 000 

0 

ICO. coo 

128 

40 

9C. 330 

23 

98 . 46 C 

2 

99. 822 

0 

100.000 

c 

ICO. ooo 

129 

52 

90.46C 

17 

98.502 

3 

99.03C 

c 

ICC. ooo 

c 

100.000 

130 

56 

9C.6CC 

27 

98.570 

2 

95. 835 

c 

ICC. 000 

c 

ICO. ooo 

131 

36 

9C.C9C 

23 

98.627 

1 

95.837 

0 

100.000 

0 

100.000 

132 

47 

9C.607 

25 

98.690 

5 

99.850 

c 

1CC.00C 

c 

ICO. ooo 

133 

40 

5C.SC7 

14 

98.725 

0 

99.850 

c 

ICO. 000 

c 

LCC.000 

134 

40 

91.C27 

26 

98.790 

6 

99.865 

0 

100.000 

0 

ICO. ooo 

135 

40 

91.127 

20 

98.040 

5 

99.877 

0 

IOC. 000 

0 

ICO. ooo 

136 

42 

91.232 

18 

98.885 

0 

99. 077 

0 

ICO. 000 

0 

ICO. ooo 

137 

51 

91.36C 

17 

98.927 

3 

99.885 

c 

ICO. 000 

c 

100.000 

138 

48 

91.48C 

14 

98.962 

3 

99.892 

0 

1CC.OOC 

c 

ICO. ooo 

139 

40 

91.58C 

15 

99.000 

4 

99.902 

0 

100.000 

c 

ICO. ooo 

140 

54 

91.715 

1 1 

99.C27 

3 

99. 910 

0 

100.000 

c 

ICO. ooo 

141 

38 

91.810 

23 

99.085 

2 

99.915 

0 

ICC. OOO 

c 

100.000 

142 

36 

91.9CC 

8 

99.105 

0 

99.915 

0 

IOC. 000 

c 

100.000 

143 

45 

92. C 12 

13 

99.137 

2 

99.92C 

0 

ICO. ooo 

0 

100.000 

144 

27 

92.C8C 

9 

99.160 

3 

99.927 

0 

1CO.OOC 

c 

ICO. ooo 

145 

40 

92. 18C 

1 1 

99.187 

C 

99.927 

0 

ICC. ooo 

c 

ICO. ooo 

146 

43 

92.287 

17 

99.230 

1 

99.930 

0 

1C0.00C 

0 

100. coo 

147 

48 

92.4C7 

7 

99.247 

1 

99.932 

0 

ICC. ooo 

c 

ICO. ooc 

148 

43 

92.515 

15 

99.285 

3 

99. 94 C 

0 

loc.ooa 

0 

100.000 

149 

33 

92.597 

12 

99.315 

3 

99.947 

0 

1CC.00C 

c 

ICC. coo 

150 

39 

92.695 

13 

99.347 

2 

99.952 

0 

ICC. 000 

c 

LOO. OOO 

151 

52 

92. 6 25 

9 

99.370 

0 

99.952 

0 

100.000 

0 

ICO. coo 

152 

37 

92.917 

10 

99.395 

0 

99.952 

0 

100.000 

c 

1C0.C0C 

153 

38 

93.C12 

10 

99.420 

4 

99.962 

0 

ICC. ooo 

c 

ICO. ooo 

154 

33 

93.C95 

5 

99.432 

C 

99.962 

0 

1CC.00C 

c 

ICO. ooo 

155 

33 

93.177 

8 

99.452 

0 

99.962 

c 

ICC. ooo 

0 

100.000 

156 

27 

93.245 

9 

99.475 

3 

99.970 

c 

100. ooc 

c 

ICC. coo 

157 

44 

93.355 

7 

99.492 

2 

99.975 

c 

ICO. ooo 

0 

ICO. ooo 

158 

38 

S3.45C 

9 

99.515 

0 

99.975 

c 

100.000 

c 

100.000 

159 

41 

93.552 

3 

99.522 

2 

99.980 

0 

100.000 

0 

ICO. ooo 

160 

40 

93.652 

6 

99.537 

0 

99.980 

0 

IOC. ooo 

c 

ICO. ooo 

161 

33 

93.735 

7 

99.555 

0 

99.980 

0 

100. ooc 

c 

10C.C00 

162 

45 

93.847 

9 

99.577 

0 

99.980 

0 

ICO. ooc 

c 

ICO. ooc 

163 

29 

93.920 

7 

99.595 

1 

99.982 

0 

ICC. 000 

c 

100.000 

164 

44 

94 . C 30 

9 

99.617 

c 

99.982 

0 

ICO. ooo 

c 

1C0.C0C 

165 

35 

94. 117 

3 

99.625 

l 

99.985 

0 

ICC. ooo 

0 

ICO. ooo 

166 

29 

94. 19C 

3 

99.632 

1 

99.987 

0 

ICC. ooo 

c 

ICO. coo 

167 

36 

94.28C 

3 

99.640 

0 

99.987 

0 

ICC. ooo 

c 

ICO. ooo 

168 

36 

94.37C 

C 

99.640 

1 

99.990 

0 

100.000 

c 

ICC. coo 

169 

49 

94.492 

9 

99.662 

0 

99.990 

0 

ICC. OOO 

0 

ICO. ooo 

170 

24 

94.552 

3 

99.670 

0 

99.990 

G 

ICC. OOO 

c 

100.000 

171 

40 

94.652 

3 

99.677 

0 

99.990 

0 

100.000 

c 

ICC. coo 

172 

26 

94.717 

5 

99.690 

0 

99.990 

0 

10C.0CC 

c 

100.000 

173 

29 

94 . 79 C 

5 

99.702 

0 

99.990 

0 

IOC. ooc 

c 

ICO. ooo 
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174 

29 

94.E62 

2 

99.707 

0 

99.990 

c 

ICC. 000 

0 

ICO. coo 

175 

31 

94.940 

3 

99.715 

0 

99.990 

0 

IOC. 000 

c 

ICC. 000 

176 \ 

30 

95.C15 

3 

99.722 

0 

99.990 

0 

100.000 

0 

100.000 

177 

30 

95.C5C 

2 

99.727 

0 

99.990 

0 

IOC. 000 

c 

ICO. 000 

178 

26 

95.155 

8 

99.747 

0 

99.990 

c 

IOC. 000 

c 

ICO. 000 

179 

31 

95.232 

3 

99.755 

0 

99.990 

0 

IOC. 000 

c 

ICO. coo 

180 

31 

95.310 

4 

99.765 

1 

99.992 

c 

ICO. 000 

0 

ICO. 000 

181 

29 

95.382 

2 

99.770 

0 

99.992 

c 

100.000 

c 

ICO. 000 

182 

26 

95.447 

2 

99.775 

0 

99.992 

0 

100. coo 

0 

ICO. 000 

183 

33 

95.53C 

4 

99.785 

0 

99.992 

0 

ICC. 000 

0 

ICO. 000 

184 

38 

95.625 

9 

99.807 

0 

99.992 

0 

ICC. 000 

c 

ICO. 000 

185 

19 

95.672 

2 

99.812 

0 

99.992 

0 

100.000 

c 

100.000 

186 

35 

95.76C 

2 

99.817 

0 

99.992 

0 

100.000 

0 

ICO. 000 

187 

29 

95.832 

1 

99.820 

0 

99.992 

0 

IOC. 000 

c 

ICC. 000 

188 

25 

95.895 

3 

99.827 

0 

99.992 

0 

ICC. 000 

c 

ICO-OCC 

189 

22 

S5.95C 

2 

99.832 

1 

99.995 

0 

ICC. 000 

c 

100.000 

190 

30 

96.C25 

3 

99.840 

1 

99.997 

0 

100.000 

0 

ICG. COO 

191 

39 

96. 122 

1 

99.842 

c 

99.997 

0 

100. ooc 

c 

100. ooc 

192 

28 

96.192 

C 

99.842 

0 

99.997 

0 

100.000 

c 

100.000 

193 

38 

96.287 

1 

99.845 

0 

99.997 

0 

ICO. 000 

0 

100.000 

194 

22 

96.342 

1 

99.847 

0 

99.997 

c 

ICO. 000 

c 

ICC. 000 

195 

24 

96.4C2 

c 

99.847 

0 

99.997 

0 

10C-000 

0 

ICO. 000 

196 

38 

96.497 

l 

99.850 

0 

99.997 

0 

100.000 

c 

ICO. 000 

197 

50 

96.622 

2 

99.855 

0 

99.997 

0 

100.000 

0 

100. ooc 

198 

35 

96.710 

1 

99.857 

0 

99.997 

0 

100.000 

0 

1CC.OOO 

199 

37 

96*80 2 

4 

99.867 

0 

99.997 

c 

100.000 

0 

ICO. 000 

200 

33 

96.885 

2 

99.872 

0 

99.997 

0 

100.000 

c 

ICO. OOC 

201 

26 

96.S5C 

1 

99.075 

0 

99.997 

0 

100.000 

c 

ICO. 000 

202 

24 

97. CIO 

3 

99.882 

0 

99.997 

0 

lOC.OOC 

c 

100.000 

203 

25 

97.C72 

2 

59.887 

0 

99.997 

0 

ICO. 000 

c 

ICC. coo 

204 

18 

97.117 

2 

59.892 

0 

99.997 

0 

100.000 

0 

100.000 

205 

21 

97.17C 

3 

59.900 

0 

99.997 

0 

100.000 

0 

100.000 

206 

21 

97.222 

1 

99.902 

0 

99.997 

0 

100.000 

0 

100.000 

207 

19 

97.27C 

C 

99.902 

0 

99.997 

0 

ICC. 000 

a 

100.000 

208 

20 

97.320 

c 

59.902 

0 

99.997 

0 

IOC. 000 

c 

ICO. 000 

209 

16 

97.36C 

2 

99.907 

0 

99.997 

0 

100.000 

c 

100.000 

210 

17 

97.402 

l 

59.910 

1 

ICC. 000 

0 

IOC. 000 

0 

100.000 

211 

1 1 

97.430 

1 

99.912 

0 

100.000 

0 

100. ooc 

c 

ICO. 000 

212 

26 

97.495 

C 

99.912 

c 

100.000 

0 

IOC. 000 

0 

ICO. 000 

213 

21 

97.547 

3 

99.920 

c 

IOC. 000 

0 

100.000 

0 

100.000 

214 

25 

97.610 

1 

59.922 

c 

100.000 

c 

ICC. ooc 

c 

ICC. 000 

215 

13 

97.642 

2 

59.527 

0 

100.000 

0 

100.000 

0 

ICO. 000 

216 

19 

5 7.69C 

0 

59.927 

0 

100.000 

0 

IOC. 000 

c 

ICO. 000 

217 

19 

97.737 

1 

99.930 

c 

100.000 

0 

IOC. 000 

c 

ICO. coo 

218 

21 

97.790 

2 

99.935 

c 

100.000 

0 

IOC. 000 

0 

ICO. 000 

219 

1 8 

97.835 

3 

59.942 

c 

100. coo 

G 

ICC. ooc 

0 

100.000 

220 

21 

97.887 

5 

99.955 

0 

100.000 

0 

IOC. 000 

0 

ICO. 000 

221 

19 

97.535 

C 

59.955 

0 

100.000 

0 

100.000 

o 

ICC. 000 

222 

22 

97 • 99 C 

l 

59.957 

0 

100.000 

0 

100.000 

0 

ICO. 000 

223 

21 

98.C42 

C 

99.957 

0 

100.000 

c 

ICO. 000 

0 

100.000 

224 

13 

98.C75 

2 

99.962 

0 

100.000 

c 

100.000 

c 

IOC. 000 

225 

14 

98. 11C 

C 

99.962 

0 

IOC. 000 

0 

ICClOOO 

c 

ICO. 000 

226 

27 

98.177 

c 

99.562 

0 

100.000 

0 

ICO. 000 

' 0 

100.000 

227 

18 

98.222 

c 

59.962 

0 

IOC. 000 

0 

ICO. 000 

c 

ICO. COO 

22 8 

16 

58.262 

L 

59.965 

0 

100.000 

0 

100.000 

c 

ICO. ooc 

229 

21 

98.315 

c 

99.965 

c 

lOC.OOC 

c 

ICC. 000 

0 

ICO. 000 

2 30 

14 

90.350 

l 

59.967 

0 

10C.000 

c 

lac.ooo 

c 

ICO. 000 

231 

9 

98.372 

1 

59.970 

0 

IOC. 000 

0 

IOC. 000 

c 

1C0.0CC 

2 32 

11 

58.4CC 

c 

99.970 

0 

100.000 

0 

ICO. 000 

0 

ICO. 000 

233 

18 

58.445 

c 

99.970 

c 

10C.000 

0 

lOC.OOC 

0 

ICO. COO 
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234 

13 

98*477 

0 

55.970 

C 

ICC.OOC 

c 

IOC. ooc 

c 

100.000 

235 

13 

98.51C 

1 

99.972 

0 

100.000 

0 

IOC. 000 

c 

100.000 

236 

18 

98.555 

1 

59.975 

0 

100.000 

0 

ICO. 000 

c 

100.000 

237 

12 

98.585 

0 

99.975 

c 

IOC. 000 

0 

IOC. 000 

0 

ICC. coo 

238 

13 

98.617 

1 

99.577 

0 

IOC. 000 

0 

100.000 

c 

ICC. coo 

239 

6 

98.632 

0 

59.977 

0 

10C.000 

0 

ICC. 000 

0 

100.000 

240 

17 

98.675 

1 

59.980 

c 

100.000 

c 

ICC. 000 

c 

ICO. 000 

241 

9 

58.697 

1 

99.982 

0 

100.000 

0 

ICO. 000 

c 

ICO. 000 

242 

5 

98.710 

0 

99.982 

0 

IOC. 000 

c 

ICC. 000 

C 

100.000 

243 

13 

98.742 

1 

99.985 

0 

100.000 

0 

IOC. OOC 

c 

ICO. 000 

244 

15 

96.78C 

1 

99.987 

0 

100.000 

0 

IOC. 000 

C 

ICO. 000 

245 

12 

98.eiC 

c 

59.987 

0 

10C.00C 

c 

100.000 

0 

100.000 

246 

9 

98.832 

c 

99.987 

0 

10C-000 

0 

ICC.OOC 

0 

100 . ooc 

247 

10 

58.857 

c 

59.987 

0 

100.000 

0 

100.000 

0 

100.000 

248 

5 

58.87C 

1 

59.99C 

0 

10C.C0C 

0 

ICC.OOC 

0 

ICO. 000 

249 

11 

98.897 

0 

59.990 

0 

100.000 

0 

ICC.OOC 

c 

ICO. 000 

2 50 

10 

58.922 

c 

59.990 

0 

IOC. 000 

0 

100.000 

0 

100.000 

251 

11 

98.95C 

c 

55.990 

c 

IOC. 000 

0 

100.000 

c 

100.000 

252 

10 

98.975 

0 

99.990 

0 

100.000 

0 

100 . ooc 

0 

ICO. 000 

253 

7 

58.952 

c 

99.990 

c 

100 . coo 

0 

ICO. 000 

c 

ICO. 000 

2 54 

12 

99.C22 

1 

59.992 

0 

100 . ooc 

0 

100.000 

c 

100.000 

255 

11 

99. CSC 

0 

59.992 

0 

100.000 

c 

ICC. 000 

0 

ICO. 000 

256 

9 

99.072 

c 

99.992 

0 

L0C.000 

0 

ICC.OOC 

c 

100.000 

257 

8 

99.C92 

c 

99.992 

c 

IOC. 000 

0 

ICC. 000 

c 

ICO. 000 

258 

6 

95.107 

C 

59.992 

0 

IOC. 000 

0 

ICC. 000 

0 

100.000 

259 

7 

99.125 

c 

59.992 

c 

IOC. 000 

0 

IOC. 000 

c 

ICO. 000 

2 60 

12 

99.155 

c 

59.992 

0 

IOC. 000 

0 

ICC. 000 

0 

100 . coo 

261 

5 

99. 167 

c 

59.992 

c 

IOC. 000 

0 

100.000 

c 

100.000 

262 

9 

59. 19C 

c 

99.992 

c 

100.000 

c 

100.000 

C 

100.000 

263 

6 

59.205 

c 

59.992 

0 

100.000 

0 

100.000 

0 

ico.coc 

264 

7 

99.222 

c 

99.992 

c 

IOC. 000 

0 

IOC. 000 

0 

100.000 

265 

9 

99.245 

c 

99.992 

0 

100.000 

0 

IOC. 000 

c 

100.000 

2 66 

5 

99.257 

c 

99.992 

0 

100.000 

0 

100.000 

c 

100 . ooc 

267 

5 

99.27C 

c 

99.992 

0 

IOC. 000 

0 

ICO. 000 

c 

100.000 

268 

3 

99.277 

c 

59.992 

0 

100.000 

0 

ICC. 000 

0 

ICO. 000 

269 

6 

59.292 

c 

99.992 

0 

IOC. 000 

0 

100.000 

0 

ICC. 000 

2 70 

4 

99.3C2 

1 

59.995 

0 

IOC. ooo 

0 

100.000 

0 

ICO. 000 

271 

6 

55.317 

0 

99.995 

0 

10C. 000 

c 

ICC.OOC 

0 

ICO. 000 

272 

7 

95.335 

c 

59.995 

0 

IOC. 000 

0 

IOC. 000 

0 

100.000 

273 

5 

99.347 

0 

59.995 

0 

IOC. 000 

0 

ICO. 000 

c 

ICO. 000 

274 

5 

55.36C 

c 

59.995 

0 

100.000 

0 

IOC. OOC 

c 

100.000 

275 

6 

59.375 

c 

99.995 

0 

IOC. 000 

0 

100.000 

c 

ICO. 000 

276 

2 

99.380 

0 

99.995 

0 

IOC. 000 

0 

ICC. 000 

0 

1 GO • 000 

277 

7 

55.397 

0 

99.995 

c 

100.000 

c 

ICC. 000 

c 

ICO. coo 

276 

6 

95.412 

1 

99.997 

0 

IOC. 000 

c 

100.000 

c 

ICO. 000 

279 

6 

59.427 

G 

99.997 

0 

IOC. 000 

0 

100.000 

0 

ICC. 000 

280 

4 

95.437 

0 

59.997 

0 

100.000 

c 

IOC. 000 

c 

ICC. 000 

281 

4 

95.447 

c 

59.997 

c 

IOC. 000 

0 

100.000 

c 

ICC. 000 

2 82 

2 

99.452 

0 

99.997 

0 

IOC. 000 

C 

100.000 

0 

100.000 

283 

6 

55.467 

c 

99.997 

0 

100.000 

0 

ICC. 000 

c 

ICO. coo 

284 

5 

99.480 

c 

59.997 

0 

IOC. 000 

0 

100.000 

c 

ICO. 000 

285 

5 

59.492 

0 

59.997 

0 

IOC. 000 

0 

ICC. 000 

0 

ICO. 000 

2 66 

3 

55.5CC 

0 

59.997 

0 

100.000 

c 

ICC. 000 

c 

ICO. coo 

287 

4 

99.510 

c 

59.997 

0 

100.000 

0 

100.000 

c 

LOO. 000 

288 

4 

99.520 

0 

99.997 

c 

IOC. ooc 

c 

IOC. 000 

0 

ICO. 000 

289 

3 

99.527 

0 

99.997 

0 

100.000 

c 

ICC. 000 

c 

ICO. 000 

290 

2 

99.532 

c 

99.997 

0 

100.000 

0 

100. OOC 

0 

100.000 

291 

2 

99.537 

c 

59.997 

0 

100.000 

c 

IOC. 000 

0 

LOO. 000 

292 

1 

59.54C 

c 

59.997 

0 

100.000 

0 

IOC. 000 

0 

ICC. 000 

293 

2 

99.545 

0 

59.997 

0 

IOC. 000 

0 

100.000 

0 

100.000 
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294 

3 

59,552 

C 

99.997 

C 

IOC. 000 

0 

100.000 

c 

ICO. coo 

295 

4 

95.562 

C 

99.997 

c 

IOC. 000 

0 

ICC. coo 

c 

ICO. 000 

296 

3 

59.57C 

c 

99.997 

0 

100.000 

0 

ICC. 000 

c 

ICC. coo 

297 

1 

99.572 

0 

99.997 

0 

100.000 

0 

100.000 

C 

lOC.OOO 

298 

7 

55.590 

c 

99.997 

c 

10C.000 

0 

100.000 

c 

ICC. coo 

299 

l 

59.592 

0 

99.997 

0 

IOC. 000 

0 

ICO. 000 

0 

ICC. 000 

300 

6 

59.607 

0 

99.997 

0 

IOC. 000 

c 

ICC. 000 

c 

ICO. 000 

301 

4 

99.617 

0 

99.997 

c 

100.000 

0 

100.000 

c 

100.000 

302 

5 

95.630 

0 

99.997 

0 

IOC. 000 

0 

lOC.OOO 

0 

100.000 

303 

6 

59.645 

c 

99.997 

0 

IOC. 000 

c 

1QC.C0C 

c 

ICO. 000 

304 

1 

99.647 

0 

99.997 

C 

IOC. 000 

0 

100.000 

c 

ICC. coo 

305 

L 

55.65C 

0 

99.997 

c 

100.000 

0 

lOC.OOO 

c 

ICO. 000 

306 

1 

95.652 

0 

59.997 

c 

100.000 

G 

100.000 

0 

ICO. 000 

307 

1 

55.655 

c 

99.997 

0 

100.000 

0 

100.000 

0 

ICO. coo 

308 

3 

99.662 

1 

ICO. 000 

0 

IOC. 000 

0 

lOC.OOO 

c 

100.000 

309 

4 

99.672 

0 

ICC. 000 

0 

IOC. 000 

c 

100.000 

c 

lOC.OOO 

310 

1 

99.675 

0 

ICO. 000 

0 

100.000 

0 

100.000 

0 

100.000 

311 

2 

99.68C 

0 

ICC. 000 

c 

IOC. 000 

0 

100.000 

0 

ICO. 000 

312 

3 

95.667 

0 

ICC. 000 

0 

100.000 

c 

ICC. 000 

c 

ICO. 000 

313 

2 

59.652 

c 

100 . coo 

0 

100.000 

0 

lOC.OOO 

0 

100.000 

314 

4 

55.7C2 

0 

ICC. coo 

0 

lOC.OOO 

0 

ICO. 000 

c 

100.000 

315 

5 

95.715 

c 

ICO. coo 

0 

100.000 

c 

ICC. 000 

c 

ICC. 000 

316 

5 

59.727 

c 

ICC. 000 

0 

lOC.OOO 

0 

ICO. 000 

C 

100.000 

317 

4 

95.737 

c 

ICC. 000 

0 

lOC.OOO 

c 

100.000 

0 

100.000 

3 18 

4 

55. 747 

0 

ICC. 000 

0 

100.000 

c 

ICO. 000 

c 

100. coc 

319 

1 

5 9* 7 5 C 

c 

ICC. 000 

0 

lOC.OOO 

0 

100.000 

0 

ICO. 000 

320 

0 

95.750 

0 

ICO. 000 

0 

100.000 

0 

1CC.00C 

c 

100.000 

321 

1 

99.752 

c 

ICC. 000 

0 

100.000 

c 

ICC. 000 

c 

ICC. coo 

322 

1 

99.755 

c 

ICC. coo 

c 

100.000 

0 

100.000 

c 

ICO. 000 

323 

4 

95.765 

0 

ICC. COO 

0 

lOC.OOO 

0 

100.000 

c 

ICO. 000 

324 

0 

59.765 

0 

ICC. 000 

0 

100.000 

0 

100.000 

c 

100.000 

325 

C 

99. 765 

c 

1 Co . 000 

0 

100.000 

0 

ICO. 000 

0 

100.000 

326 

2 

99.77C 

0 

ICC. 000 

C 

lOC.OOO 

0 

lOC.OOO 

c 

ICO. 000 

327 

0 

99.770 

c 

ICC. 000 

0 

100.000 

c 

ICC. 000 

c 

lOC.OOO 

328 

0 

99.77C 

c 

ICC. 000 

c 

lOC.OOO 

0 

lOC.OOO 

0 

100.000 

329 

0 

'99.770 

c 

ICO. 000 

c 

lOC.OOO 

0 

100.000 

c 

ICO. 000 

3 30 

0 

99.770 

0 

100. COO 

0 

LOC.OOO 

0 

ICO. 000 

c 

100.000 

331 

2 

99.775 

c 

ICO. 000 

0 

lOC.OOO 

0 

100.000 

0 

ICO. 000 

3 32 

' 1 

99. 777 

0 

ICO. 000 

0 

100.000 

0 

100.000 

0 

100.000 

333 

0 

99.777 

c 

ICC. 000 

0 

lOC.OOO 

0 

lOC.OOO 

c 

ICO. 000 

334 

2 

99.782 

c 

ICC. 000 

0 

lOC.OOO 

0 

ICO. 000 

0 

100.000 

335 

1 

99.705 

c 

ICO. COO 

0 

lOC.OOO 

0 

ICO. 000 

0 

ICO. 000 

336 

1 

99.787 

0 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ICO. 000 

337 

0 

99.767 

c 

ICC. 000 

0 

100.000 

0 

1C0.00C 

C 

100.000 

338 

2 

59.792 

0 

ICO. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. coo 

339 

3 

99.0CC 

0 

ICO. coo 

0 

lOC.OOO 

0 

ICC. 000 

0 

ICC. 000 

340 

2 

99.8C5 

c 

ICO. 000 

0 

100.000 

0 

ICC. 000 

c 

100.000 

341 

1 

59.8C7 

c 

ICC. 000 

0 

100.000 

0 

100.00C 

0 

ICO. 000 

342 

2 

99.812 

0 

ICO. 000 

0 

lOC.OOO 

0 

100.000 

C 

100.000 

343 

1 

99.815 

c 

ICC. 000 

0 

100.000 

0 

ICO. 000 

0 

100.000 

344 

1 

99.817 

0 

ICC. 000 

0 

lOC.OOO 

0 

lOC.OOO 

0 

ICO. 000 

345 

1 

99.82C 

c 

IOC. 000 

0 

lOC.OOO 

0 

ICC. 000 

0 

ICO. coo 

346 

1 

99.822 

c 

ICC. 000 

0 

100.000 

0 

lOC.OOO 

C 

ICC. coo 

347 

0 

99.822 

0 

100.000 

0 

lOC.OOO 

0 

100.000 

c 

100.000 

348 

2 

99.827 

0 

ICC. 000 

c 

lOC.OOO 

0 

lOC.OOO 

c 

100.000 

349 

0 

99.827 

0 

ICC. 000 

0 

100.000 

0 

ICO. 000 

0 

ICO. 000 

3 50 

2 

99.832 

c 

ICC. 000 

0 

100.000 

0 

1CC.000 

0 

100.000 

351 

2 

99.837 

c 

ICC. 000 

c 

100.000 

0 

100 . ooc 

0 

ICO. 000 

3 52 

0 

99.837 

c 

ICC. 000 

0 

lOC.OOO 

0 

100.000 

0 

ICO. 000 

353 

2 

99.842 

c 

ICC.OCO 

0 

lOC.OOO 

0 

100.000 

0 

ICO. 000 


TABLE VI. - Continued. COMPUTER OUTPUT OF SAMPLE PROBLEM 


354 

1 

99 • £45 

C 

ICO. 000 

0 

IOC. 000 

c 

100.000 

0 

100.000 

35 5 

3 

99*852 

0 

ICO. 000 

0 

100.000 

c 

100.000 

c 

ICO. 000 

356 

1 

99.655 

0 

ICC. 000 

0 

100.000 

0 

IOC. 000 

c 

ICC. 000 

357 

1 

99.857 

c 

ICC. 000 

0 

IOC. 000 

0 

ICO. 000 

0 

100.000 

358 

3 

99. 665 

0 

ICO. 000 

0 

IOC. 000 

0 

ICC. 000 

0 

100.000 

359 

0 

99.865 

0 

ICO. 000 

0 

100.000 

0 

IOC. 000 

c 

100.000 

360 

0 

99.665 

c 

ICC. 000 

0 

100.000 

0 

100.000 

0 

lOC.OOC 

361 

1 

99.667 

0 

ICO. 000 

0 

100.000 

0 

ICC. 000 

0 

ICO. 000 

362 

3 

99. e75 

0 

ICC. coo 

0 

100.000 

c 

ICC. 000 

c 

ICC. coo 

363 

1 

99.877 

c 

ICC. 000 

0 

IOC. 000 

0 

100.000 

c 

100.000 

364 

0 

99.877 

c 

ICO. 000 

0 

IOC. 000 

0 

ICO. 000 

c 

100.000 

365 

1 

99.880 

0 

ICO. 000 

0 

100.000 

0 

lOC.OOC 

c 

ICO. OOC 

366 

3 

99.887 

c 

ICC. 000 

0 

IOC. 000 

0 

100.000 

0 

100.000 

367 

1 

99. 690 

0 

ICO. 000 

c 

100.000 

c 

ICC. 000 

c 

100.000 

368 

2 

99.695 

c 

ICO. 000 

c 

IOC. 000 

0 

ICC. 000 

c 

ICC. 000 

369 

1 

99.897 

c 

1CC.00C 

c 

IOC. 000 

0 

100.000 

0 

ioo.ooo 

370 

l 

99.9CC 

0 

ICC. 000 

0 

IOC. 000 

0 

100.000 

c 

100.000 

371 

2 

99.9C5 

c 

ICO. coo 

0 

IOC. 000 

c 

ICC. 000 

c 

100.000 

372 

0 

99.9C5 

c 

ICC. 000 

0 

IOC. 000 

0 

ICC. 000 

c 

ICO. 000 

3 73 

l 

99.907 

0 

ICC. 000 

c 

100.000 

0 

ICC. 000 

0 

ICO. 000 

374 

1 

99.910 

c 

ICC. 000 

0 

100.000 

0 

1CC.00C 

c 

ICC. 000 

375 

0 

99.910 

c 

ICC. 000 

0 

IOC. 000 

0 

IOC. 000 

c 

ICO. 000 

376 

l 

99.912 

0 

ICC. 000 

0 

IOC. 000 

c 

ICO. 000 

c 

ICO. 000 

377 

3 

99. 92 C 

0 

1CC.0C0 

0 

IOC. 000 

0 

100.000 

0 

ICO. 000 

3 78 

C 

99.92C 

c 

ICO. 000 

0 

100.000 

0 

100.000 

0 

ICO. 000 

379 

l 

99.922 

c 

ICC. coo 

0 

IOC. 000 

0 

ICC. 000 

0 

100.000 

380 

0 

99.922 

c 

ICC. 000 

0 

IOC. 000 

0 

ICC. 000 

c 

ICO. coo 

381 

0 

99.922 

0 

ICC. 000 

0 

100.000 

0 

lOC.OOC 

c 

100 . oco 

382 

0 

99.922 

c 

ICC. 000 

0 

100.000 

0 

100.000 

0 

ICO. 000 

383 

1 

99.925 

c 

100.000 

c 

IOC. 000 

c 

IOC. 000 

c 

ICO. 000 

384 

2 

99.93C 

0 

ICC. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. coo 

385 

2 

99.935 

c 

ICC. 000 

0 

IOC. 000 

c 

ICO. 000 

c 

ICO. coo 

386 

l 

99.937 

0 

ICC. 000 

0 

100.000 

0 

ICC. 000 

c 

1C0.CC0 

387 

0 

99.937 

c 

ICC. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. coo 

388 

2 

99.942 

0 

ICC. 000 

0 

100.000 

0 

IOC. 000 

0 

ICO. 000 

389 

C 

99.942 

0 

ICO. 000 

0 

10C.000 

0 

ICC. 000 

c 

ICO. 000 

390 

1 

99.945 

c 

ICC. 000 

0 

100.000 

c 

IOC. 000 

c 

100.000 

391 

0 

99.945 

c 

ICC. 000 

0 

IOC. 000 

c 

ICC. 000 

0 

100.000 

392 

1 

99.947 

0 

ICC. 000 

c 

100.000 

0 

ICC. 000 

0 

ICO. 000 

393 

0 

99.947 

c 

ICO. 000 

c 

IOC. 000 

0 

100. OOC 

C 

100.000 

3 94 

0 

99.947 

0 

ICO. 000 

0 

IOC. coo 

c 

IOC. 000 

c 

ICO. oco 

395 

1 

99.950 

c 

ICC. 000 

0 

IOC. 000 

0 

lOC.OOC 

c 

ICO. 000 

396 

0 

99.95C 

c 

ICC. 000 

0 

IOC. 000 

0 

IOC. 000 

0 

100.000 

397 

0 

99.950 

0 

ICC. coo 

0 

100.000 

0 

ICC. OOC 

c 

ICC. coo 

398 

0 

99.95C 

c 

ICC. coo 

0 

IOC. 000 

0 

ICC. 000 

0 

ICC. OOC 

399 

1 

99.952 

c 

ICO. 000 

0 

IOC. 000 

0 

100.000 

0 

IOC. oco 

400 

0 

99.952 

0 

ICC. 000 

c 

100.000 

0 

100.000 

c 

ICC. 000 

401 

1 

99.955 

c 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ICO. 000 

402 

c 

99.955 

c 

ICC. 000 

0 

L0C.000 

0 

ICC. OOC 

c 

ICO. 000 

403 

0 

99.955 

0 

ICO. 000 

c 

IOC. 000 

0 

lOC.OOC 

c 

100.000 

404 

0 

99.955 

c 

ICC. coo 

0 

100.000 

0 

IOC. 000 

c 

ICO. 000 

405 

0 

99.955 

c 

ICC. 000 

0 

100.000 

0 

IOC. 000 

0 

100.000 

406 

0 

99.955 

c 

ICC. 000 

0 

100. 000 

0 

IOC. 000 

0 

ICO. 000 

407 

0 

99.955 

c 

ICC. 000 

c 

100.000 

c 

ICC. OOC 

c 

100.000 

408 

2 

99. 96 C 

0 

ICO. 000 

0 

100.000 

0 

lOC.OOC 

0 

ICO. 000 

409 

1 

99.962 

c 

ICC. 000 

0 

100.000 

0 

100.000 

0 

ICO. 000 

4 10 

0 

99.962 

c 

ICC. 000 

c 

IOC. 000 

0 

lOC.OOC 

0 

100.000 

411 

0 

99.962 

0 

ICC. 000 

0 

100.000 

c 

ICO. 000 

c 

100.000 

412 

c 

99.962 

c 

ICO. 000 

0 

IOC. coc 

0 

ICO. 000 

c 

100.000 

413 

2 

99.967 

0 

ICO. 000 

0 

100.000 

0 

IOC. 000 

c 

100.000 
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TABLE VI. - Continued. COMPUTER OUTPUT OF SAMPLE PROBLEM 


414 

0 

99.967 

C 

ICO. 000 

0 

10C.00C 

0 

100.000 

C 

ICO. 000 

415 

0 

99.967 

0 

ICC. 000 

0 

100.000 

0 

ICO. 000 

C 

ICO. 000 

416 

0 

99.967 

0 

ICC. 000 

0 

IOC. 000 

0 

100.000 

0 

100.000 

417 

0 

99.967 

0 

ICC. 000 

c 

100.000 

0 

ICO. 000 

c 

ICO. 000 

418 

0 

99.967 

0 

ICO. coo 

c 

100.000 

0 

ICO. 000 

C 

ICO. 000 

4 19 

0 

99.967 

c 

ICC. 000 

0 

100.000 

0 

100.000 

0 

100.000 

420 

0 

99.967 

0 

ICO. 000 

c 

IOC. 000 

0 

100.000 

0 

ICO. 000 

4 21 

0 

99.967 

0 

IOC. 000 

0 

lOC.OOO 

0 

ICO. 000 

0 

100. COO 

422 

0 

99.967 

0 

ICC. 000 

0 

100.000 

0 

lOC.OOO 

c 

100.000 

423 

1 

99.97C 

c 

ICC. 000 

0 

lOC.OOO 

c 

lOC.OOO 

0 

ICO. 000 

424 

0 

9S.97C 

0 

ICO. 000 

0 

100.000 

0 

lOC.OOO 

C 

100. COO 

425 

0 

99.S7C 

0 

ICC. 000 

0 

lOC.OOO 

0 

100.000 

0 

ICO. coo 

426 

l 

99.972 

0 

ICO. 000 

0 

lOC.OOO 

0 

100.000 

0 

100.000 

427 

0 

99.972 

0 

ICC. 000 

C 

100.000 

0 

100-000 

0 

ICO. coo 

428 

c 

99.972 

c 

ICC. 000 

0 

100.000 

0 

100.000 

0 

100.000 

4 29 

0 

99.972 

0 

ICC. 000 

0 

100.000 

c 

lOC.OOO 

0 

100.000 

430 

0 

99.972 

0 

ICC. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. coo 

431 

0 

99.972 

c 

ICC. 000 

0 

lOC.OOO 

0 

100.000 

c 

ICO. 000 

4 32 

c 

99.972 

0 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ICO.OOC 

433 

1 

99.975 

c 

ICO. 000 

0 

lOC.OOO 

0 

lOC.OOO 

0 

100.000 

4 34 

1 

99.977 

0 

ICO. 000 

0 

lOC.OOO 

0 

lOC.OOO 

G 

ICC. 000 

435 

0 

99.977 

0 

ICO. 000 

0 

100.000 

0 

lOC.OOO 

c 

100.000 

436 

c 

99.977 

0 

ICC. coo 

0 

100.000 

0 

100.000 

c 

lOC.OOO 

437 

0 

99.977 

0 

ICC. 000 

0 

lOC.OOO 

0 

100.000 

c 

100.000 

438 

c 

99.977 

0 

ICC. 000 

C 

ICO. 000 

0 

ICO.OOC 

0 

100. coc 

439 

L 

99.980 

0 

ICO. 000 

0 

100.000 

0 

lOC.OOO 

c 

100.000 

440 

0 

99.980 

0 

ICC. 000 

0 

lOC.OOO 

0 

lOC.OOO 

0 

ICO. 000 

441 

0 

99.980 

0 

ICO. 000 

0 

lOC.OOO 

0 

l'CC.000 

0 

ICO. 000 

442 

0 

99.90C 

0 

ICC. COO 

0 

100.000 

0 

100.000 

0 

ICO. 000 

443 

1 

99.982 

0 

100.000 

0 

lOC.OOO 

0 

lOC.OOO 

c 

100.000 

444 

0 

99.982 

c 

ICO. 000 

0 

100.000 

0 

100.000 

c 

100.000 

445 

0 

99.982 

G 

ICC. 000 

0 

100.000 

0 

lOC.OOO 

C 

ICO. 000 

446 

0 

99.982 

0 

ICC. 000 

0 

lOC.OOO 

0 

ICO. coo 

C 

100. oco 

447 

0 

99.982 

0 

ICC. 000 

0 

100.000 

0 

100.000 

C 

100.000 

448 

1 

99.985 

c 

ICC. 000 

0 

lOC.OOO 

c 

lOC.OOO 

C 

ICC. 000 

449 

c 

99.985 

0 

ICC. 000 

0 

100.000 

0 

lOC.OOO 

0 

ICO.OOC 

4 50 

c 

99.985 

0 

100.000 

0 

lOC.OOO 

Q 

10C.0C0 

C 

ICO. 000 

451 

0 

99.985 

c 

ICC. 000 

c 

lOC.OOO 

0 

ICC. 000 

c 

ICO. 000 

452 

0 

99.985 

c 

ICC. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. 000 

453 

0 

99.985 

c 

ICC. 000 

0 

100.000 

0 

100. 000 

c 

ICC. 000 

454 

0 

99.985 

c 

ICO. 000 

c 

lOC.OOO 

c 

100. ooc 

c 

100.000 

455 

c 

99.985 

0 

ICO. 000 

0 

lOC.OOO 

0 

lOC.OOO 

c 

100.000 

456 

0 

99.985 

c 

ICC. 000 

0 

100.000 

0 

100.000 

c 

100.000 

457 

0 

99.905 

0 

ICC. coo 

0 

100.000 

C 

ICC. 000 

c 

100.000 

458 

0 

99.905 

0 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ICO. 000 

459 

0 

99.985 

c 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ICO.OOC 

460 

0 

99.985 

0 

100.000 

0 

IOC. coo 

0 

100.000 

0 

ICO. 000 

461 

0 

99.985 

0 

ICC. coo 

0 

lOC.OOO 

C 

100.000 

0 

ICO. 000 

462 

0 

99.985 

0 

ICC. 000 

0 

100.000 

0 

100.000 

C 

ICO. 000 

463 

1 

99.987 

c 

ICC. coo 

0 

lOC.OOO 

c 

IOC. ooc 

0 

100.000 

464 

0 

99.987 

c 

ICO. 000 

0 

100.000 

0 

lOC.OOO 

0 

ICO. 000 

465 

0 

99.987 

0 

ICC. coo 

0 

lOC.OOO 

0 

100.000 

C 

ICC. coo 

466 

0 

99.987 

0 

1CC.OOO 

0 

lOC.OOO 

0 

100.000 

c 

ICO. coo 

467 

1 

99.990 

0 

ICO. 000 

0 

100.000 

0 

100.000 

0 

ICO. 000 

468 

0 

99. 990 

0 

ICO. 000 

0 

lOC.OOO 

0 

lOC.OOO 

G 

100.000 

469 

0 

99.990 

0 

ICC. 000 

0 

100.000 

0 

ICC. 000 

0 

ICO. 000 

470 

0 

S9.99C 

0 

ICO. COO 

0 

100.000 

0 

100.000 

c 

ICO. 000 

471 

0 

99.990 

0 

ICO. 000 

0 

lOC.OOO 

0 

ICO. 000 

C 

ICO. 000 

472 

c 

99.990 

0 

ICO. 000 

0 

lOC.OOO 

0 

100.000 

c 

ICO. coo 

4 73 

0 

99.990 

c 

ICO. 000 

0 

lOC.OOO 

0 

ICO. 000 

c 

100 . coo 
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474 

1 

99*992 

C 

ICC. 000 

0 

IOC. 000 

c 

ICC. 000 

c 

icc.coc 

475 

C 

59.992 

0 

ICC. 000 

0 

ICC. COO 

0 

ICC. 000 

0 

ICO. 000 

476 

c 

59.592 

0 

ICO. 000 

c 

ICC. 000 

0 

ICC. 000 

c 

ico.ooc 

477 

0 

55.592 

C 

ICC. 000 

0 

IOC. 000 

0 

IOC. 000 

c 

ICO. coo 

478 

0 

59.552 

C 

1CC.00C 

0 

IOC. 000 

0 

100.000 

c 

ICO. 000 

479 

0 

59.992 

C 

ICC. 000 

0 

100.000 

0 

1CC.C0C 

c 

ICO. COO 

480 

0 

59.552 

0 

ICC. 000 

0 

100.000 

0 

ico.ooc 

C 

ICC.COC 

48 1 

c 

59.952 

0 

1CC.C0C 

0 

100.000 

0 

100.000 

C 

100.000 

482 

c 

59.552 

c 

ICC. 000 

0 

IOC. 000 

0 

1CC.00C 

C 

100.000 

483 

c 

59.992 

c 

ICC. 000 

0 

100.000 

0 

100.000 

c 

ico.ooc 

484 

0 

59.592 

c 

ICC. 000 

0 

IOC. 000 

c 

KC.OOC 

c 

ico.ooc 

485 

0 

59.592 

0 

ICC. 000 

0 

IOC. 000 

0 

10C.C0C 

c 

ICO. coo 

486 

i 

59.995 

c 

ICO. 000 

0 

100.000 

0 

ICC. 000 

c 

ICO. 000 

487 

i 

59.597 

c 

ICC. 000 

0 

IOC. 000 

0 

100.000 

c 

ICO. 000 

488 

0 

59.997 

c 

100.000 

c 

IOC. 000 

0 

ICO. 000 

c 

ICC.COC 

489 

0 

59.597 

0 

ICC. coo 

0 

IOC. 000 

0 

ICC. 000 

c 

ICC. 000 

490 

0 

59.997 

0 

ICC. 000 

c 

IOC. 000 

0 

ICO. 000 

c 

100.000 

491 

0 

59.997 

0 

1CC.00C 

0 

100.000 

0 

100. ooc 

c 

ICO. coo 

492 

0 

55.597 

c 

ICC. COO 

0 

100.000 

0 

100.000 

c 

ICO. coo 

493 

0 

59.997 

c 

ICC. coo 

0 

100.000 

c 

100.000 

c 

ICO.OOC 

494 

0 

99.997 

0 

ICC. 000 

0 

100.000 

C 

ICO. 000 

0 

ICC.COC 

495 

0 

59.997 

0 

ICC. 000 

0 

IOC. 000 

c 

100.000 

0 

ICO. COO 

496 

0 

59.597 

0 

ICO. 000 

c 

IOC. 000 

C 

ICC. 000 

0 

ICO. 000 

497 

0 

59.557 

c 

ICC. 000 

0 

100.000 

c 

ICC. 000 

c 

ICC. 000 

498 

c 

59.597 

c 

ICC. coo 

0 

IOC. 000 

c 

ICC. 000 

c 

ICC. 000 

499 

c 

99.997 

c 

ICC. COO 

c 

100.000 

0 

IOC. 000 

0 

100.000 

500 

c 

59.997 

c 

ICC. 000 

0 

IOC. 000 

c 

ICC. 000 

c 

ICO. 000 

INflNlTY 

1 

1CC.CC0 

c 

ICC. coo 

0 

IOC. 000 

c 

300.000 

c 

ICO. 000 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 24, 1967, 
124-09-01-03-22. 
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APPENDIX A 


RADIUS AND DISTANCE EFFECT CALCULATION AT ANY TIME DURING MISSION 

In this calculation, the mission is assumed to be composed of a transfer ellipse 
from the Earth to a planet, a hold time at the planet, and a transfer ellipse to return 
from the planet to the Earth. Both transfer ellipses have the Sun at one focus. One 
such transfer ellipse is shown in figure 3. 



The mission requires a transfer from the Earth at point A to the planet at point B 
along the ellipse. Some parameters of the ellipse are as follows: Points E and F are 
foci, D is the perihelion point, DG is the major axis, the eccentricity e = EF/ DG, the 
semilatus rectum P = CE (which is perpendicular to the major axis), also 
P = (1 - e 1 2 )DG/2. 0, and the true anomaly 0 is the polar angle measured from the line 
between the Sun and perihelion. 

The polar representation of the ellipse (with the Sun as the origin) is 


1 + e cos 0 

which relates the distance r from the Sun and the true anomaly 0 for known values of 
e and P. 

The time elapsed in going from the perihelion position (0 = 0) to any 0 is (e. g. , 
see ref. 5) 
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where 

H central force constant, GM 
M mass of central body 

G universal gravitational constant = 6.668x10 km /(sec )(g) 

For heliocentric orbits: 

M = Mq = Mass of Sun = 1. 989X10 00 g 

M = GMq = 1. 32X10 11 km 3 /sec 2 = 9. 906xl0 20 km 3 /day 2 

The time t(0 Q ) corresponding to the true anomaly at the start of the mission is de- 
termined, and the elapsed time to reach any position in the transfer is 

T(e) - t(e) - t(e D ) 

To relate the distance r to an elapsed time T requires determining the true 
anomaly 9 corresponding to T by using equation (A2), and then the distance r corre- 
sponding to 6 by using equation (Al). Because 6 cannot be obtained explicitly from t 
in equation (A2), a numerical solution is necessary which requires a disproportionate 
amount of computation time. To improve computational efficiency, the code evaluates 
r corresponding to elapsed time T by the following procedure. 

The total angle traversed on the transfer ellipse, 9 Q to 6^ is divided into 49 equal 
intervals so that 


A 6 


e f~ 


49 


and 


9^ = 8 Q + i A9 

The values of T(0 j) are determined for each of the 50 angles directly from equation (A2) 
and the distances determined for each of the angles directly from equation (Al). A sim- 
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ilar set of values of 8 , T, and r are obtained for the second transfer ellipse to return 
from the planet to the Earth. Then, interpolating between these values of T and r, a 
table of the distance effect g(r) = (l/r) a at the beginning of each day of the trip is con- 
structed. To find g(r) at any time T during the trip, the computer determines T to 
the nearest day and a g(r) corresponding to that day simply by noting its value at this 
position in the table. The g(r) for the first day is the first one in the table, the g(r) 
for the second day is the second entry in the table, and so on. Inasmuch as r varies 
slowly with time, this method of selecting g(r) does not introduce appreciable error. 
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APPENDIX B 

UNCERTAINTY IN DETERMINATION OF D p 

For a given trip, let the actual dose probability density be f(D) so that f(D)dD rep- 
resents the true probability that the dosage for a trip is in the interval dD about D. 

The cumulative probability that a trip has a dosage less than D is given by 

rD 

p(D)= / f(?)d? (Bl) 

J o 

Equation (Bl) yields p as a function of D (see fig. 4) 

One is concerned with the inverse problem, that is, to find the dose corresponding 
to a given p. Let the true value of D corresponding to a given p be denoted by D*, 
as shown in figure 4. 




Figure 4. - Frequency and cumulative probability distributions. 
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It is desired to determine how much an estimate Dp (of D*) might be in error when 
the estimate is based on a group of M identical trips simulated by the computer code. 

The estimate D is taken to be the dosage encountered on the (pM)** 1 trip, where 
the trips have first been arranged in order of increasing dosage (it is assumed that the 
number pM is an integer). Thus, in the group of M trips, (pM - 1) of the trips en- 
counter doses less than Dp and (1 - p)M = qM of the trips encounter doses larger than 
D . The probability that the measured value of D lies in the interval dD lying at D 
is given by the trinomial expression 

f D (D )dD = M! ( p + 5) pM_1 (q - 6) qM f(D JdD (B2) 

U p P (pM - 1)1 (qM)l p 


where 


P + 



(B3) 


is the true probability that the dosage encountered on a trip is less than Dp. 

The situation is as depicted in figure 4. Note that equation (B3) may be broken down 
into the two equations: 


P = 



(B3a) 


6 = 



(B3b) 


If it is assumed that 5 is much smaller in magnitude than either p or q and that 
numbers (pM - 1) and qM are large with respect to p or q, equation (B2) may be 
written in approximate form as 


Vv 


,Mf 2 < Vexp/-i M «‘ 


27rpq 


2 pq 


(B4) 


If D is sufficiently close to D*, equation (B3b) may be approximated by 
P P 
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« = f(D p )(D p - Dj) 

which when substituted into equation (B4) yields 




) = 



(B5) 


The distribution, as given by equation (B5), is a normal distribution with a mean value 
equal to D* and a standard deviation equal to 

ir 



i 


pq 


Mf 2 (Dp) 


t 


P(1 - P) 

M^ 2 (Dp) 


(B6) 


Some of the algebra involved in arriving at equation (B4) from equation (B2) follows: 
From equation (B2) 

In fp. (DJ = (pM - l)ln(p + 6) + qM ln(q - 6) + In Ml - ln(pM - 1)1 - ln(qM)l + In f(D ) 
u n P F 


This equation can be expanded as follows: 


(pM - l)ln(p + 5) = (pM - l)ln p 



= (pM - l)ln p + (pM - 


= pM In p - In p + (pM 
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qM ln(q - 6) = qM In q 


= -qM In 


H) 

— + qM In ( 1 - — \ 

q V V 


-qM In — - qM 

q 


« + i(»\ 2 +i/«Y + 

q 2\q / 3\q/ 


Using Stirling's formula for large M, 


Ml 


1 s V2 h 


m M+(V2) 


M 


one may express 


Similarly 


and 


In Ml = — In 2tt + ( M + ~ )ln M - M 

2 V 2/ 


= - In 2tt + M In M + - In M - M 


-ln(pM - 1)1 = - - In 2n - ( pM - - )ln(pM - 1) + (pM - 1) 

2 V 2/ 


-ln(qM)l = 


- In 273- - (c 

2 \ 


qM + - ) ln(qM) + qM 
2 / 


= - — In 277 - qM In q - qM In M - -i In q - — In M + qM 


Adding yields 
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In Ml - ln(pM - 1)1 - ln(qM)l = - — In 2ir + pM In M - pM ln(pM - 1) 

2 


+ * ln(pM - 1) + qM In - + - In - - 1 
2 q 2 q 


Also ln(pM - 1) can be expanded as follows: 
ln(pM - 1) = ln(pM) | 


'(l ) 

\ ml 


= In p + In M - | + - — + — - + 


V pM 2(pM) 2 3(pM) 3 


Combining these terms results in 


In 


f D (°d) = - - ^ 2tt + pM[ J- + 1 + L_ 

p 2 \ pM 2(pM) 2 3(pM)‘ 


+ • • 


+ — lnp + — InM - -(— + + + . 

2 2 


2 \ pM 2(pM) 2 3(pM) 3 


+ - In - - l-lnp + (pM - 1) 
2 q 


p 2 \p / 3 \p / 


- qM 


« + iM 2 + iM 3 + 

q 2\q / 3\q / 


+ In 


£(Dp) 
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and combining some terms 



Further simplifying yields 
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Assume that 6 is an order of magnitude less than either p or q and M is sufficiently 
large so that all terms in the brackets other than 


pM/s\ 2 _ qM /g\ 2 
2 \p/~ 2 \q/ 


can be neglected. Then 



1 i n _2vrpq_ _ 

2 Mf 2 (D p ) 


which leads to equation (B4). 


1 M5 2 

2 pq 
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APPENDIX C 


STRATIFIED SAMPLING OF FIRST FLARE OCCURRENCE DURING TRIP 

To reduce the uncertainty in Dp for a given number of trips, the occurrence of the 
first flare is forced to follow the exponential distribution by a stratified sampling tech- 
nique (see ref. 6). This is accomplished as follows: The random number interval (0, 1) 
is subdivided into M equal parts (where M is the number of trips evaluated in a group). 
The midpoint of the first interval = 1/2M is the random number used to determine 
(from eq. (1)) when the first flare occurs oh the first trip, the midpoint of the second in- 
terval £ 0 = 3/2M is the random number used to determine when the first flare occurs 

Ct 4-U 

on the second trip, and, in general, = (2i - 1)/2M for the i trip. When a ^ is ob- 
tained so that the time of occurrence exceeds the total trip time, no flares occur on that 
trip. After the time-to-first-flare is selected, the procedure for selecting the flare 
type and time to subsequent flares is as indicated in the METHOD OF ANALYSIS section. 
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APPENDIX D 


FORTRAN IV LISTING 

The complete IBM 7094-n FORTRAN IV listing of the MCFLARE program is pre- 
sented here. The random number generator, RAND/ SAND (DECK RANGPL) included 
here will execute properly only on a 36-bit word machine such as the IBM 7090/7094/ 
7094-11 with a FORTRAN IV compiler. If this program is used on a different computer, 
the user will have to supply his own random number generator and alter the two CALL 
RAND statements and one CALL SAND statements to suit. 


SIBFTC RANGPL DECK 

SUBROUTINE SAND ( IR ) 

C A RANDOM NUMBER GENERATOR FOR A 36 BIT MACHINE WITH 

C MAXIMUM INTEGER CAPABILITY OF 2**35-l 

I R = 1 

ZR= 1. 0/(2. 0**35 -1.0) 

K= 5**15 
RETURN 

ENTRY RAND ( R ) 

I R = I R*K 

R= FLOAT ( IR) * ZR 

RETURN 

END 
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SIO *YuD1520 G. P. LAHTI— MCFLARE 
$1 EJOB 

S I B FTC MCFLAR CECK 

OUSE CALCULATION FCR SCLAR FLARES 
G.P.LAFTI NA S A— LE Vt I S RESEARCH CENTER 

MONTE CARLC FLARE OCCLRRENCE CALCULATION 
ALL LOSE BINS ARE UNIT BINS 

EACH EVENT IS ASSURED TO OCCLR CNLV ONCE DURING THE TIRE TPEF 
MAX- NUMBER OF FLARF TYPES = 25 

MAX. NUMBER OF SHIELDS = 7 

MAX. NUMBER OF DOSE BINS = 9CC 

COMMON NTYPE. NT* NBIN* NB, NSHELD * TOTALT* UMTS* PP* PPNT* 

1 PZFREC(2), POFREO ( 1C* 2 ) * IDI36)* SH I ELD ! 7 ) 

DIMENSION I FREGZ ( 26 ) * I FREC ( 50 *26 ) * JBIN(7,901) 

DIMENSION J FREOZ ! 26 ) * JFREC ( 50 .26 ) * KBIN(7,901) 

DIMENSION CDOSE 125*7) * NUMI26), TCTOOSI7)* TCURI25J* R(IOOC) 

500 FURMAT16I5.2E1C.4) 

501 FORMAT! 8FS. 4) 

502 FUKMATt 12A6 ) 

503 FORMAT! 1FC/ 1HC/19HCT0TAL TRIP TIME = F8.2.5H CAYS) 

50 A FORMAT! 1 H 1 ✓ ! 1H .12A6)) 

505 FORMAT! 1 HC . 1 6 * 4 6H * C CC TRIPS TC BE RUN HITH FRINTCUT AFTER EACH, 

115, 4 H , C C C J 

506 FORMAT !29FC DOSES BASED CN l./R**,I2) 

663 FURMAT ! 28HL THIS CALCULATION C C N S I DER S , I 3 . 1 2H FLARES, AND,I3*16H S 
IHIELLIS. NAMELY* 7I2X.A6)/ 42H0 OOSE CALCULATION USES 

2LJP TO. 1 5 , 1 OH UNIT BINS ) 

£64 FGRMAT118FL FLARE INPUT DATA, 25X* 9F0CSES IN *A6// 

167H FLARE OCCURRENCES OF DURATION COSE FROM ITH TYPE F 

ALARE / 

267H TYPE. ITH TYPE FLARE CF ITH TYPE AT 1 AU MTH JTH SHI 

HELD- / 

34 X * IHI *5X,2HINF7. 1* 19H DAYS FLARE* DAYS, 7!6X,A6)) 

665 FORMAT! 15* J8X, IHI. E14.4* 7E12.4) 

695 FOR MA T ( 0 J 2 ) 

CALL SANCIXYZ) 

READ! 5,659 ) XYZ 
PUNCH 695, XY l 

1 READ! 5,5C; J ! ID! I )* 1=1,36) 

READ! 5.5CC ) MTHOl S . NTHOUS. I XP ,N T YPE . NSHELD ,NB I N , TREF,TCTALT 
MTHUUS T HC l SAND TRIPS ARE RUN 

INTERMEDIATE PRINTOUT IS MADE AFTER EACH NTHOUS THOUSAND TRIPS 
DISTANCE EFFECT (FROM SUN) IS ASSUMEO TO BE 1/R**IXP 
IF! ( NTFCUS.LE.C) .OR. ( NTHCU S .G T . M THCU S ) ) NT H0US= MTH OUS 

MT IMFS = MTI-OUS/NTHCUS 
M THOU S = MT IMES*NTHOLS 

NTYPE = NUMBER CF DIFFERENT TYPES CF FUARE EVENTS CQNSICEREC 
NSH FL 0= NUMBER CF DIFFERENT SHIELDS CONSIDERED 
NBIN = NUMBER OF DOSE BINS 

TREE = LENGTH OF TIME OF REFERENCE PERICD FOR INPUT FLARE CCCURRENCES 
T 0 1 AL T = TOTAL TRIP TIME* DAYS (USED ONLY IF IXP=0> 

RE A 0 <5*5C;j UNITS, (SHIELDCI1* 1 = 1,7) 

UNITS = COSE UNITS USED (E.G. *RADS*) (USED IN TITLING OUTPUT) 

SH I EL D= SHIELD DE S CR I P T I C N 1 E CR OUTPUT TITLING ONLY) 

WRITE (6. 5C4 l( ID! I ). I = 1*36) 


° n on no 


WR I TE ( 6 « 6 63 ) NTYPE*NSHELD,(SHIELD(I) *1 = 1*7) • NBIN 
WRI Tt(6, 5C5 )MTHOUS,NTHCUS 
WRITb(6,5C6) IXP 

WRITE 16,6641 UM1S, TREE, (SHIELD! I i • 1=1,7) 

DO 3 I = 1 • NT YPE 

3 READ (5. 50) !DOOSE(I,J), J=1.NSHELD) 

C DOOSFU.JJ = DOSE AT 1 AU FROM ITH TYPE FLARE, WITH JTH SHIELC 

READ! 5*50 J (TDLR(I), I=i,NTYPE) 

C TDJR!I) = TIME DURATION OF ITH TYPE EVENT 

DO 4 I=1*NTYPE 

4 WRITE (6,665) I, TDUR(I), (DDOSE(I.J), J=l,NSHELD) 

C START CALCLLATION 

IF! IXP >6,5,6 

5 WRITF(6,5C3) TOTAL T 

IF! TOTALT.IE. C.C) CALL EXIT 
JDAYS=TOTALT+L.C 
DO 105 1= 1,JDAYS 
105 R ( 11*1.0 
GO TO 7 

6 CALL RADIISIIXP, TOTALT, R) 

7 NT = NTYPE+I 
QNT=NTYPE 

C ON T = TO TAL NO OF EVENTS WHICH OCCUR IN TIME TREF 

PP= TOTALT/TREF 

C PP = FX PEC T E C OCCURRENCES OF EACH FLARE TYPE CN TRIP 

PPNT=PP*CNT 

C PPNT=EXPECTED TOTAL NUMBER OF FLARE EVENTS TO CCCUR CN TRIP 

ISTAR = TREF/ONT 

POISSON THEORETICAL PROBABILITIES FOR C TC 10 OCCURRENCES 
PER TRIP OF EACH FLARE TYPE 
PROB=PP/CNT 
DO 13 J= 1 , 2 
PR0B=PRG8*CNT 
PZFRFG(J) = EXP(-PROB) 

PZEREO! J )=PROBABIL ITY OF ZERO FLARES OCCURRING 
PU FR EO ( 1 , J ) = PZFREQ (J )*PRCB 
00 13 1=2, 1C 
F I = I 

FREQUENCIES ARE FOR I OCCURRENCES OF THE JTH TYPE FLARE 
13 PUFRE Q ( I , J ) = POFREO! I—1»J)*PRCB/FI 
NB = NBIM1 
00 16 1=1, NB 
00 16 J = 1 , N SHELD 

16 KBIN! J, I ) = 0 
DO 17 1=1, NT 
JFRFGZ ( I ) =C 
DO 17 J = 1 , 5 C 

17 JFREQI J, I) =0 
KLJMP = 1CCC+NTH0US 
XT=2*KLUMF 
XT=1.0/XT 
XXT=XT+XT 

MAIN ITERATION LOOP 
00 99 NTIMES =1 .MTIMES 
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C 


18 


19 


00 18 I = 1 , NB 
DO 18 J= 1 » N SHELD 
JB I N ( Jf I ) = C 

JBINI J, I) = NO. OF TRIPS 
00 19 1=1, NT 
IFRFQZtl) =c 
no 19 j=i.«c 
I FR EG ( J , I ) = C 
IFR EGZ ( I ) = NO. GF TRIPS 
IFREOI J, I ) =N0. OF TRIPS 
NOTE — 1= N T Y PE+ 1 = NT 
I FREU ( J * NT ) = NO. OF TRIPS 


IN ITH DOSE BIN WITH JTF SHIELD 


IN WHICH NO I TYPE FLARE OCCURREC 
IN WHICH ITH TYPE FLARE CCCURREC JT IMES 
FEANS ANY TYPE FLARE 

WITH J FLARES OCCURRING (ALL TYPES) 


XN=-XT 


INTER MECIATE LOOP 


DO 60 INDEXN=1.NTH0US 
DO 60 I NCE XT= 1 » 1000 
00 21 I = 1 , N T 

21 NUM I I ) = C 

C NUM ( I ) = NC. OF TIMES ITH TYPE FLARE OCCURS IN TRIP 

C NUM(NT) = TOTAL NC. OF FLARES THAT OCCURRED ON TRIP 

DO 22 1= 1, NSFELD 

22 r0TD0S(I)= C. 

C TOT DO SI I )= DOSE FOR ITH SHIELD 

DAYS = 0. 

XN- XN ■fX X T 

C SET UP FIRST FLARE FCR TRIP FROM STPATIFIEC SAMPLE 

PR AND = XN 
GO TO 232 

23 CALL R ANC ( PR AND ) 

233 T = TSTAR*(-ALOG(PRAND) ) 

DAYS = DAYS+T 

IF! DAYS.GT. TOTALT) GO TO 31 
NUMINT) = N UM ( NT ) + 1 

C SELECT FLARE TYPE OCCURRING IN THIS INTERVAL 

CALL RANCIPRAND) 

C FOR EACH FLARE OCCURRING ONLY CNCE IN THE REFERENCE PERIOD 

I=GNT*PRAND-H.O 

C AN I TYPE FLARE HAS OCCURRED 

26 NUM I I ) =NUM I I ) + 1 

C CALCULATE DISTANCE EFFECT FACTOR • RR 

JDAY=DAYS* 1 .0 
RR=RI J C A Y ) 

C ITH TYPE FLARE HAS OCCURRED — MAKE DOSE CALC FCR ALL SHIELDS 

DO 28 K=1.NSHELD 

28 TOrDOS(K>= TOTDOS(K) + CDOSE(I.K) *PR 
DAYS = DAYS + TDURI I ) 

I F ( DAYS— TCTALT ) 23,31,31 

C 31 END OF TRIP ANALYZE RESULTS 

31 I FI NUM (NT) ) 3 3 . 3 3 , 41 

C 33 NO FLARES HAVE OCCURRED ON THIS TRIP 

33 DO 34 1=1. NT 

34 IFREOZI I) = IFREG Z ( I ) + 1 
DO 36 J = 1 .NSHELD 

35 JBINI J, 1 ) = JBINt J , 1 ) +1 
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GO TO 60 

C 41 SONF FLARES HA VE OCCURRED GN THIS TRIP 
4 1 DO 50 1 = 1 *NT 

I F 1 NUN I I ) i 4 2, 42* 43 

42 IFREQ7 ( I ) = IFREOZ (I) +1 

C NO FLARES CF TYPE I HAVE OCCURRED CN THIS TRIP 

GO TO 50 

43 NN = NUR ( I ) 

C NN FLARES OF TYPE I HAVE OCCURRED ON THIS TRIP 

I FR SO ( NN * I ) = I FREC I NN * I )+l 
50 CONTINUE 

DO 54 J=1*N$HEUD 
PUT TOTAL COSES IN CORRECT BINS 
FOR UNIT CCSE BINS..*. 

I = TOTDOS{JJ-H.O 
IFU-GT.NBJ I =N8 
54 J B I N ( J* I )= JBINU»I I +1 
60 CONTINUE 

CALL FINALE! IFRECZ* IFRFQ, JBIN* KLUMP) 

00 6 5 I = 1 * N 8 
DO 65 J=J*N$HELD 

65 KBIN( J* I ) = KB I N ( J * I ) 4- JBIN(J*I) 

DO 70 1 = 1 *NT 

JFRFOZII) = JFREOZ (11+ IFRECZ(I) 

DO 70 J= 1 * 5 C 

70 JFRFO(J.I) = JFREQ<J,I) + IFREC(J*I) 

I F ( NTIMES . EQ • 1) GO TO 99 
KUUMPS = N ^ INES*KUUMF 

CAUL FINALEt JFRECZ. JFRE 0 * K0IN, KLUHPS) 

99 CONTINUE 

PUNCH 699* XYZ 
GO TO 1 
END 


tIBFTC FINAL 


CECK 


SUBROUTINE F ! NAL E ( I FR ECZ , I FREC * JBIN* MANY) 

COMMON NT YPE , NT, NB I N • NB. NSHELD, TOTALT. UNITS, PP. PPNT, 

1 PZFREQ12). POFREQ! 1C, 2) • IDI36), SH I ELD ( 7 ) 

DIMENSION I FREQ Z < 26 ) , I FRE C ( 50 ,26 ) , JBIN(7,901) 

DIMENSION K SUM I 26 ) ♦ CPCT(7.9C1> 

600 FORMAT! I 3, 1618) 

601 FORMAT! 5FCSUMS, 16, 1518) 

602 FORMAT (112, AX .7116. F8. 3)) 

603 FORMAT i AX , 8H I NF I N I T Y , AX, 7(16, F8, 3)) 

660 FORMAT ( 6 1 M THE NUMBER OF TRIPS IN NHICH I FLARES OF TYPE J OC 

1C UR — / 6 6 F (THE LAST COLUMN IS THE NUMBER OF TRIPS IN WFI 

2CH I FLARES OF ANY TYPE OCCUR)/ 1 IHC I J = i , 1518) 

661 FORMAT ( 2CF 1 DOSE DISTRIBUTION /1H0 , 15X ,7 ( 8H ShIELC-A6 

1)// 1 6X « 7 ! 1 AH NC. PERCENT)/ 5X,5HD0SE , A6, 7(1AH IN LESSTH 

2AN ) / 1 6h UPPER LIMIT, 7!1AH BIN LIMIT )) 

662 FORMA T ( A 2FC THE EXPECTED NC. CF EACH TYPE IS ...F8.3/ 

1 A 2F THE EXPECTED NC. GF ANY TYPE IS ...F8.3) 

671 FORMAT ( 73HL CALCULATED PROBABLE OCCURRENCES PEP TRIP USING POIS 

ISON DISTRIBUTION /35HC NUMBER EACH TYPE ANY TYPE/ 

29H0 C. 1P2E13.5) 

672 FORMAT! 19, 1P2E 13.5 ) 

DO 6A I = 1 , 5 C 

IMAX = 51-1 

IF ( IFREO! IMAX, NT) )6A. 6A. 66 
6A CONTINUE 
66 DO 68 J = 1 , NT 

KSUM ( J ) = I FREOZ ( J ) 

DO 68 1 = 1, IMAX 

68 KSUM! J ) = K S LM ( J ) + IFPEC(I.J) 

C SEARCH DC S E BINS FOR HIGHEST OCCURRENCE 

Nbd = NB+1 
INFINY =1 
DO 70 1=1, NB 
MAX B= NBB-I 
DO 70 J = 1 , N SHE LD 
IF! JBIN(J.MAXB) ) 7C.7C.72 
70 CONTINUE 

72 1 F ( I.EG.n INF I N Y = 2 
FMISS = MANY 
SS= 100. /FMISS 
DO 74 J = 1 . N SHELD 
NS JM = C 

00 7 A I = 1 , MAXB 

NSUM = NSUM + JBIMJ.I) 

C SUM= N SUM 

7 A CPCT! J • I )= CSUM+SS 
NM= MINO! 16. NT) 

WRITF(6,66C)( J, J=2.NM) 

1 = 0 

HR I TE (6.6CC) I.( IFFECZ(J) ,J=1*NM) 



on 76 I* I. IMAX 

76 WRITE ( 6 • 6 C G ) I • CIFPEfl 1 1 * J ) * J*1 .NM ) 

WRITE ( 6 « 6C 1 ) ( K SUM (J) , J = 1 , NM ) 

I F I NT • L E • 1 6 3 GO TO 77C 
MR I TE ( 6 * 663 ) < J * J=18,NT) 

663 FORMA T( 1 1 EC I J = 17, 1518) 

1 = 0 

WRITE ( 6 • 6 C C ) 1,1 IFPEC7 (J) *J=17,NT) 

00 7661=1* I RAX 

766 WRITE { 6 , 6CG ) I * I IF PEC (I , J ) , J= 1 7 , NT ) 

WRITE ( 6 * 6C 1 ) (KSUM(J) *J=17.NT) 

770 WR I TF { 6 » 662 ) PP . PPM 
C OUTPUT THEORETICAL OCCURRENCES PER TRIP 

WRITE(6.671) PZFREC < 1 ) , PZFREG(2> 

WRITE <6*672)1 I, POFREGII.il. PCFREC(I*2). 1=1*10) 
WRITE16.661) (SHIELD(I), 1=1,7). UNITS 
GO TO (78,77), INFINV 

77 MX= MAXB -1 
GO TO 8 C 

78 MX = MA XB 

80 00 82 1=1, MX 

82 WRITE (6,602) I, (JBIN(J.I). CFCTIJ.I). J=1,NSHELD) 

GO TO (99*84) » INF INY 

8 4 WR I TE ( 6 . 6C 3 ) ( JB IN ( J ,M AXB ) , C PCT ( J , MAX 8 ) , J=l *NSHEL D) 

99 RETURN 
END 
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SIBFTC RAOMCF 


CECK 


SUflKOUTINF PAD ILS( I XP t N • TOTAL!, RRP) 

DIMENSION T HET A ( 2 ) , THETAF ( 2 ) , E(2), P(2). AA(2), TAU(2>, OAU(IOC) 
1 , riME(lCC). W(l), R R { IOC) , PRRUOOO) 

DATA P I • A L , U/ 3. 1415*527. 1.4959E8, 9.90598E20/ 

U = PRODUCT OF MASS OF SUN AND UNIVERSAL GRAVITATION CONSTANT 
U HAS UNITS OF KM**3/DAY**2 

AU = MEAN CISTANCE BETWEEN EARTH AND SUN. KM 
601 FORMAT! F 6 . C , 2 { 3F 6 . C . E 9 . 4 ) ) 

75C FORMA T ( 2 CM MISSION PARAMETERS/ 

123H0 TOTAL TRIP TIME. DAYS. F9.4// 

21 7X.28FEARTF TO PLANET, TC EARTH/ 

3 L 9 H SEMI-LATUS RFCTUM, E13.5, E13.5, 3X. 2HK M / 

414H ECCENTRICITY. 5X, 2F13.5/ 

523H TRUE ANOMALY AT S T AR T , F 9 . 3 ,F 1 3 . 3 , 3X , 7HCEGRE ES/ 

62 1 H TRUE ANOMALY AT END, F 1 1 . 3 , F 1 3 . 3 . 3 X . 7 F C E G R E F S / 

714H PERICC. DAYS.5X.2FI3.4/ 

8 1 7 H SEMI- MAJOR AXIS,E15.5,F13.5,3X,2HKM/ 

947HL ELAPSEC TIME (DAYS) VS DISTANCE FROM SUN (AU)) 

751 FORMAT! 5 ( I 5 . F 9 . 3 . F 7 . 3 ) ) 

RADIUS GENERATICN SECTION 

RF AD (5,60) THCLD, (TFETA(I). THETAF(I), E(I), P(I), 1=1,2) 

THETA = TRLE ANOMALY — ANGLE FROM PERIHELION AT START — FOR ELLIPSE 

T HE TA F= TRLE ANOMALY — ANGLE FROM PERIHELION AT END FOR ELLIPSE 

*** THETA, THETAF UNITS ARE DEGREES**** THETAF MUST RF .CT. THETA 
F = ECCENTRICITY {■= 0 FOR CIRCULAR ORBIT) FGR ELLIPSE 

P = SEMI-LATUS PECTUM IN KILOMETERS PCR ELLIPSE 

ELLIPSE 1 IS FOR EARTH TO PLANET 
ELLIPSE 2 IS FOR PLANET TC EARTH 
r I M F ( 1 )= c . 

DO 60 J = 1 , 2 
EM = 1 . — F ( J ) 

F P = 1 • * E ( J ) 

FF = F M* E P 

R = SORT (P(J)**3/L) / FE 

C = 2./SCRT(FE) 

D = SCPTIEM/EP) 

A A i J ^ =P I J ) /EE 

T A U ( J ) =2 •*FI*SCRT (AA (J)**3/L ) 

C AA = SEMT-MAJOR AXIS OF ELLIPSE 

C TAU = PERIOD OF ELLIPSE (CAVS) 

DTE E I A = (THETAF(J)— THETA(J) ) /49. *PI/180. 

X = T H E T A ( J ) *PI/IEC. 

GO TO (15,10, J 
1G T IME ( 51 )= 7 I ME ( 5 C ) + THOLC 

15 DO 50 1 = 1, 5C 

K= I J — 1 ) * 5 C +1 

DFiM = l.+E( J )* COS(X) 

I F ( A R S ( X— PI) -l.CF-B) 16,16,17 

16 7 = — P I / 2 . 

GU TO IB 



17 L = ATAMC*SIN( X/2. ) / C0SIX/2.J ) 

18 TT = B* I C*Z-E< J )*SIN< > I/DEM 

C TT = TIME ELAPSED IN TRANSFER CCNIC AS RECKONED FRCM PERIHELION 

IF! I. EC.) ) GO TO 22 
D T = T T-1TT 

THEU)= TIME(K-l) + DT 

C TIMF(K) = TIME ELAPSEC SINCE START CF TRIP 

19 IH TIMEIK) - TIME(K-l) ) 20, 2C, 22 

C DT IS SOMETIMES LCU BY A PERIOD OR ThO.. THIS CP CORRECTS THIS 

20 T I M F ( K ) = TIMEIK) + TAU<J) 

GO TO 1 9 

22 T T T = T T 

X= X+DTFETA 

DAUIK) = P ( J ) /DEN / A U 

RR ( K ) = DAL (K )**{— IXPGN ) 

50 CONTINUE 
6C CCNTINLE 

TOTAL T = 7 I ME I ICC ) 

r!RITF(6,75C ) TIME! IOC I* P(l), F(2), E(l), E { 2 ) * THETA(l), 
ITHfcTAm, THE TAF ( 1 ) t ThETAF(2), T AL I 1 ) * TAUI2), AA(1), AA(2) 

DO 70 I = 1 » 2 C 
12= I + 2 C 
13= I+«C 
14= 1+60 
15= I + 8 C 

70 WR1TE16,751) It TIME(I), DAU(I)t I2» TIMEU2), DAUII2I, 13, 

1 T I MF ( 1 3 I * DAUTI3), 14, TIME(I4), DAUII4). 15, TIME(I5), CAUII5) 
C 

C SFT UP l/CAl**IXPCN = RRRU) FCR EACH CAY IN THE TRIP 

JT IMF = TO TAL T 
N M i N = 2 

DU 99 1= 1. JTIME 

T= I 

DO 93 N=NMIN,100 

IFi T-TIME (N I i 94,53,53 

93 CQiMTJNLE 

94 N M IN = N 

R = I r — T IMEfN — 11 J*( RR(M- RR I N- U ) / I T I ME I N )- T I M E ( N- 1 ) ) + RR(N-l) 

KRRI I )=R 
99 CGNTINLE 

RRR< JTIME+1 ) = RRR( JTIME) 

RETURN 

EiNO 
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